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Introduction 


A MECA workshop, "Dust on Mars III,” was sponsored by NASA 
through the Lunar and Plane'ary Institute and held on September 
21-23, 1988. The workshop was hosted by the University of Colo- 
rado's Laboratory for Atmospheric and Space Physics at the Stanley 
Hotel in Estes Park, Colorado; logistical and administrative support 
was provided by the LPI Projects Office. Two related MEC A work- 
shops, “Dust on Mars” and “Dust on Mars II,” were held at Arizona 
State University in 1985 and 1986. 

The goal of this workshop was to stimulate cooperative research 
on, and discussion of, dust-related processes on Mars, and to provide 
valuable background informa r ion and help in planning of the Mars 
Observer mission. The workshop was organized to address the fol- 
lowing general questions: (1) How is dust ejected from the martian 
surface into the atmosphere? (2) How does the global atmospheric 
circulation affect the redistribution of dust on Mars? (3) Are there 
sources and sinks of dust on Mars? If so, where are they and how do 
they vary with time? (4) How many components of dust a)e there 
on Mars, and what are their properties? 

The program committee consisted of P. Christensen (Arizona 
State University), B. Jakosky (University of Colorado), S. Lee (Uni- 
versity of Colorado), and R. Zurek (Jet Propulsion Laboratory). To 
address the questions posed above, the committee structured the 
workshop into four primary discussion sessions: (1) dust in the 
atmosphere, (2) dust on the surface, (3) dust properties, and 
(4) dust observations from future spacecraft missions. 

Thirty people attended the meeting and participated in the dis- 
cussions. Of these, 23 made oral presentations of their dust-related 
research. This technical report provides a summary of the presenta- 
tions and discussions, and a compilation of the contributed 
abstracts. 


Steven Lee 
Workshop Organizer 
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Program 


Wednesday Morning, September 21, 1988 
8:30 a.m. 

SESSION I: DUST IN THE ATMOSPHERE: DUST STORMS, ATMOSPHERIC CIRCULATION MODELS, 

AND SEDIMENT TRANSPORT 
Chairman: R. W. Zurek 

Martian Great Dust Storms: Aperiodic Ph nomena? 

R. W . Zurek arid R . M. Haberle 

Simulations of the General Circulation of the Martian Atmosphere II. Dust Storms 
J. B. Pollack R M. Haberle, J. Schaeffei , and H. Lee 

Mars Global Atmospheric Oscillations: Annually Synchronized, Transient Normal Mode Oscillations and the Triggering of 

Global Dust Storms 
J. E. Tillman 

Numerical Simulations of Global Dust Stcim Decay 

J. Murphy, O. B. Toon, J. B. Pollack, aid R. M. Haberle 

Simulations of the Martian Boundary Lay»*r : Factors Controlling the Behavior of the Surface Stress 
R. M. Haberle,]. B. Pollack, and). Sl:Iu u’ffer 

Use of Threshold Velocities in Dust Prediction 
D. Gillette 

BarocliniC'Radiative Instability and Marti m Summer Polar Cyclones 
]. R. Barnes 

Earthbased Monitoring of the Martian Atmosphere Dust Opacity 
T. Z Martin 

Correlations Between Surface Albedo Features and Global Duststorms on Mars 
L J. Martin and P. B . James 

Comparison of Martian Aeolian Features and Results from the Global Circulation Model 
R. Greeley, A S/cypeck, and J. B. PoRaiT 

Polar Cap Edge Heating in the Dusty Non hern Winter on Mars 
R. Kahn 

Discussion 
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Thursday Morning, September 22, 1988 
8:30 a.m. 

SESSION II: DUST ON THE SURFACE: SEDIMENT DEPOSITS AND SEDIMENT TRANSPORT 

Chairman: P. R. Christensen 

Dust in the Polar Regions of Mars: Is It the Same as at Low Latitudes? 

P. C. Thomas 

The North Polar Sand Seas: Preliminary Estimates of Sediment Volume 
N. Lancaster and R. Greeley 

Martian Sediment Transport: Evidence from High Resolution Thermal Data 
P. R. Christensen 

Viking Photometric and Albedo Studies of Regional Dust Transport on Mars 
S. W. Ijee and R. T. Chancy 

^Possible Intense Vortices and the Potential for Dust and Sand Transport on Mars 
J. A. Grant and P. H. Schultz 

Discussion 


Thursday Afternoon, September 22, 1988 
2:00 p.m. 

SESSION III: PROPERTIES OF MARS DUST 
Chairman: B. M. Jakosky 

An Assessment of the Meteoric and Meteoritic Contributions to Dust on Mars 
G. J. Flynn and D. S. McKay 

The Mineral Components of Dust on Mars 
A Bunin 

Depolarized Radar Echos and Dust on Mars 
T. W. Thompson and H. J. Moore 

Geographical Variability of Mars Surface Physical Properties 
B. M. Jakosky and J . A. Marshall 

Preliminary Results and Overview of Spectral Remote Sensing Observations of Mars During the 1988 Opposition 
J. F. Bell III, T. B. McCord, P. G . L ucey, and D. L Blaney 

^Infrared Transmission Measurements of Martian Soil Analogs 
T. L Roush 


Discussion 
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Friday Morning, September 23, 1988 
8:30 a.m. 

SESSION IV: DUST OBSERVATIONS FROM FUTURE SPACECRAFT MISSIONS 

Chairman: L. W. Esposito 

Mars Observer PMiRR: Distribution and Transport of Atmospheric Dust 
R. W Zurek, R. D. Haskins, and D. }. McCIeese 

*The Mars Observer Thermal Emission S lectrometer: Studies of the Composition of Martian Dust 
P. R. Christensen 

Mars Observer VIMS and Phobos KRFM ISM: Observations of Mars Dust 
L W. Esposito 

GENERAL DISCUSSION 
WORKSHOP ADJOURNS 


*Abstracts not in abstract volume. 
^Presented by title only. 
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Workshop Summary 


DUST IN THE ATMOSPHERE: DUST STORMS, 
ATMOSPHERIC CIRCULATION MODELS, AND 
SEDIMENT TRANSPORT 

The workshop’s first session focused on research related 
to the Mars atmosphere. Discussions included results of 
modeling (such as global circulation models), observations 
of atmospheric properties and circulation patterns, sediment 
transport processes and patterns, and observations and 
phenomena related to the production and timing of dust 
storms. Eleven oral presentations were given. 

R. Zurek and R. Haberle considered the possibility that 
planetary-scale dust storms may l>e a relatively recent 
phenomenon (Martin; Icarus , 1984). Considerable debate 
ensued as to whether observations before 1950 had 
sufficient coverage to establish that global dust storms did 
not occur during that period. Tht observational records, 
both historical and modern, do show that great dust storms 
have not occurred every Mars year. In particular, the 
presentation by L. Martin and P. James showed that there 
was definitely no great dust storm m 1986. Thus the great 
dust storms cannot be counted on to close the seasonal 
cycles, as proposed for the water eye le, for instance (Davies; 
Icarus , 1981). Furthermore, a scale factor (1/b) is needed 
to go from the net annual effects of great dust storms 
in years of their occurrence to much longer time scales, 
like an obliquity cycle. At present the data indicate that 
1 < b < 10 (4-6 major dust storms <n 50 Mars years). Zurek 
and Haberle suggested that this scab* factor might be refined 
using a model that treats great dust storms as a chaotic 
dynamical system, in a manner sim lar to recent treatment 
of terrestrial interannual climatic e\ cuts like El Nino (Vallis; 
Science, 1986). 

J. Pollack, R. Haberle, J. Schaeffer, and H. Lee presented 
results from a series of steady state simulations generated 
using the NASA Ames Mars General Circulation Model 
(MGCM). These simulations were for different seasonal 
dates and for different global dust amounts, fixed in time 
and uniformly distributed in space. The presentation 
focused on two aspects of the simulations. First, maps of 
the maximum surface stress computed by the MGCM 
indicated that the wintertime, high northern midlatitudes 
and the subtropics of both hemispheres were the latitudinal 
zones of maximum surface stress. Zurek suggested that it 
would be useful to identify the diurnally varying 


contribution to the maximum surface stress, since 
atmospheric tidal models predict that planetary-scale, daily 
varying near-surface winds would be largest in the 
subtropical regions. The second point was that the presence 
of dust in the polar night increases significantly the 
atmospheric condensation of CO 2 at the expense of surface 
condensation. This atmospheric condensation is due to 
radiative cooling by the airborne dust and provides a 
mechanism for scavenging the airborne dust from the polar 
atmosphere. 

R. Greeley, A. Skypeck, and J. Pollack compared MGCM 
maps of maximum wind stress at the surface with the 
orientation and location of bright and dark streaks. 
Generally, the bright streaks correlated well with the 
MGCM simulations, except in Tharsis. The dark streaks 
did not correlate well; their distribution appears less 
organized, suggesting that the dark streaks are controlled 
more by local topographic effects. There is relatively good 
correlation of rock abundance with the MGCM predictions 
of the longitudinally varying pattern of surface wind stress, 
but not with areas of apparently very long-term erosion, 
such as lulls, which are possibly wind-sculptured. 

J. Tillman showed interannual differences in planetary- 
scale dust storms as revealed in the Viking Lander pressure 
data record. He showed that large-scale dust storms 
occurred in 1977 and in 1982, but probably not in the 
two intervening martian years. The appearance of nearly 
diurnal and semidiurnal modes during late northern summer 
(L s ~-150°) suggests the appearance of resonant atmospheric 
normal modes. The nature of much larger amplitude modes 
of nearly diurnal frequency during the onset of the 1977b 
storm is more uncertain; the issue of whether such modes 
are a cause, in conjunction with other atmospheric 
circulation components, or an effect of the dust storm onset 
is unresolved. 

J. Murphy, O. B. Toon, J. Pollack, and R. M. Haberle 
presented 1 -D simulations of dust clearing at lower 
latitudes. They found that disc-shaped particles (i.e., 
particles with an aspect ratio of 1/10) were needed to 
explain the observed rates of dust clearing, following the 
1977a,b dust storms. Spherical particles fell out much too 
quickly, even in the presence of very large (107 cm 2 /s) 
vertical mixing. The two 1977 great dust storms had 
different clearing rates that were not easily explained by 
changing the particle size distribution or the vertical eddy 
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diffusivity. The difference between the two rates may be 
due to different large-scale vertical motions at the latitudes 
of the Viking Landers during the two storms. This possibility 
is currently being investigated using a 2-D version of the 
aerosol model. 

As noted by Haberle, Pollack, and Schaeffer, there are 
better working hypotheses for how great dust storms 
originate, than for why the great dust storms almost 
immediately begin to decay after their explosive growth. 
They outlined two possibilities. First, the storms may 
exhaust the supply of dust in their source regions. This 
suggestion led to an interesting discussion of whether mobile 
dust can be sequestered locally or whether such dust can 
be literally removed from relatively large source regions. 
If the latter can occur, how is dust resupplied to these 
regions, which historically have been the growth centers 
of both local and great dust storms? Also investigated was 
the possibility that there was a local feedback mechanism 
in which the increased static stability of a dusty near-surface 
atmosphere reduced the surface wind stress, thereby 
shutting off the supply of dust to the atmosphere. Detailed 
modeling of this feedback in the planetary boundary layer 
at the subsolar point, presumably in the source region of 
the great dust storm, shows too modest a weakening of 
the surface wind stress to account for the observed time 
scales of dust storm decay. 

D. Gillette revisited the issue of determining the 
threshold friction velocity for movement of material at 
the surface. The apparent or effective velocity required 
for such movement may be 2-3 times greater than the 
computed threshold, if there are nonerodible elements 
distributed about the surface. On Earth, a knowledge of 
the size distributions of aggregates and the physical state 
of the soil can be used to compute the expected wind 
erosion by convolving the velocities actually required for 
movement on a given surface and the likelihood of those 
velocities being exceeded there (i.e., the surface wind 
climatology). It is not obvious to what degree the empirical 
relations derived from terrestrial data can be adapted to 
Mars. 

J. Barnes reexamined a mechanism proposed to explain 
the rare spiral clouds observed in summer at high northern 
latitudes. This mechanism consists of a radiative instability 
in which the vertical scale height of the dust is quite small, 
so that an initial perturbation resulting in upward vertical 
velocity will bring up more dust, leading to more solar 
heating in the dustier air and thereby producing a stronger 
vertical velocity. Using a beta-plane dynamical model and 
a linearized transport equation, Barnes investigated systems 


that were both radiatively and baroclinically unstable. He 
found that both mechanisms would greatly amplify the 
initial perturbation with e-folding growth times of just 
1-2 days. However, there was little scale selection in the 
growth rates, meaning that the observed scale of the spiral 
clouds did not grow faster in the model than did 
perturbations characterized by other horizontal scales. 

R. Kahn noted an apparent dichotomy of the zonal- 
mean wind directions at high northern and southern 
latitudes during their respective spring seasons. Wind 
directions inferred from Viking visual observations of lee 
waves indicate that westerlies continued to prevail 
equatorward of the retreating polar cap edge during early 
spring in the north (10° < L s < 40°), while there was a 
shift to easterlies equatorward of the springtime south 
seasonal cap. Kahn argued that this may be due to thinner 
than expected frost deposits equatorward of 70°N, citing 
orbiter observations of patchiness in the early spring frost 
cover and the different rates of retreat of the two seasonal 
caps. He further argued that this patchiness was due to 
greater poleward heat transport inhibiting CO 2 
condensation in northern winter and that this enhanced 
transport was due to the greater dustiness of the atmosphere 
at that time. Haberle noted that an inviscid, thermally 
direct circulation would continue to produce easterlies even 
without a cap sublimation flow and suggested that the 
observed lee waves were associated with westerly baroclinic 
waves. Kahn countered that viscous models (notably that 
of Magalhaes; Icarus, 1987) could produce zonal-mean 
westerlies. The ensuing discussion pursued the issue of 
whether or not the zonal-mean circulation at these times 
and latitudes is viscously driven, and if so, what the physical 
nature of that internal viscosity might be. 

T. Z. Martin discussed the Earth-based monitoring of 
martian atmospheric dust opacity by imaging Mars at five 
thermal IR wavelengths in the Earth’s 10-5 micron 
“window.” These images spectrally span the silicate dust 
feature used to estimate atmospheric dust optical depths 
from Mariner 6 and 7 IRS r Mariner 9 IRIS, and 
Viking IRTM observations of Mars. This technique has the 
potential to survey Mars even during those seasons when 
daytime viewing is required. Earth-based data have been 
acquired in 1986-1988, but no funds are available for the 
reduction and analysis of the images. 

As noted earlier, there was no great dust storm in 1986, 
nor had a planetary-scale storm occurred yet this year 
(L s < 270 °). Martin and James noted that in 1986 even 
local dust storm activity seemed less than in previous years. 
Dark albedo features in 1986 were generally very prominent 
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and sharp. Solis Lacus in particular was large, dark, and 
well defined; in fact, it appeared darker than at any time 
since the 1926 opposition. 

DUST ON THE SURFACE: SEDIMENT DEPOSITS 
AND SEDIMENT TRANSPORT 

The purpose of the second session was to focus on the 
physical properties of martian surface sevlimentary materials 
(ranging from dust to sand), and the transport of such 
materials. Four oral presentations addiessed these topics, 
using a variety of remote sensing techniques including 
thermal, visual, and photometric observations. 

The initial paper of the session, by P. Thomas, dealt 
with observations in Acidalia, Amaze mis, and Argyre of 
“linear” dust streaks that have remarkably sharp boundaries 
and great lengths. These streaks exceed several hundred 
kilometers in length, with a slight increase in width along 
their length, and often have a slight curvature. They are 
also remarkable in that they aie regenerated in 
approximately the same location and configuration 
following global dust storm deposition. Among the questions 
raised by these streaks are the mechanisms by which their 
sharp boundaries are maintained end whether these 
mechanisms are related to atmospheric dynamics or to 
surface properties. 

N. Lancaster and R. Greeley described detailed estimates 
of sand volumes in the martian north polar erg. The 
thickness of these deposits was determined by measuring 
the spacing of dunes and applying a teirestrial relationship 
between spacing and height. Once the 1 leight was estimated, 
an average, or equivalent sand thickness (EST), can be 
determined. The major assumptions in this estimate are 
(1) that martian dunes have similar shapes to those on 
the Earth, and (2) that the sand estimate is only for the 
sand contained in the observed bed forms. With these 
assumptions the sand estimate for the north polar erg ranges 
from an average sand thickness of 0.1 to 6.0 m, with a 
mean of 1.8 m. These estimates arc lower than those 
previously obtained, and suggest that the erg could have 
been derived from erosion of the layered deposits. To 
provide the derived amount of sand would require the 
erosion of only — 275 m of layered retrains that contained 
1% sand. 

P. Christensen outlined the cunent questions and 
possible avenues for solution of the martian sedimentary 
cycle. Of particular interest are ( 1 ) the |>ossibly cyclic nature 
of dust deposition and removal; (2) the processes that 
modify, indurate, and eventually erode dust deposits; (3) the 


period(s) of deposition and erosion; (4) how do dust deposits 
feedback int- > the general circulation; and (5) what has been 
the transport history of both dust and sand. One data set 
that may shed some new light on the questions of cyclic 
deposition and erosion is a moderate (~30 km) resolution 
thermal inettia image of Mars. This image reveals discrete 
units of differing thermal inertia in Arabia and Lunae 
Planum that are roughly concentric to current dust deposits. 
Based on tht ir location and the discrete increases in thermal 
inertia at their boundaries, they may be the remnants of 
earlier, indurated dust deposits. 

S. Lee and T. Clancy described their photometric analysis 
of Viking IRTM broadband reflectance measurements taken 
at multiple phase angles over a specified surface target during 
a single orbit. Using these observations it is possible to 
fit a photometric function, in this case using Hapke’s theory, 
to the data to derive parameters related to surface 
properties. The primary parameter that was adjusted to 
provide the best fit to the data was Hapke’s compaction 
parameter. The areas studied to date were in Acidalia 
Planitia, Amazonis, Arabia, Cerberus, and Olympus Mons. 
These preliminary results suggest that there are differences 
in compaction parameter between bright and dark surfaces, 
and between high and low inertia surfaces. In addition, 
there appear to be significant differences in this parameter 
within small regions. Future work will focus on detailed 
atmospheric corrections and the mapping of photometric 
properties over larger regions. 

PROPERTIES OF MARS DUST 

The third session concentrated on the chemical and 
physical pr iperties of martian dust, including discussion 
of laboratory and remote sensing results. Five oral 
presentations were included. 

G. J. Flynn presented a discussion of the amount of 
meteoric and meteoritic material that would be intermixed 
with surface material of martian origin. This is of interest 
in analysis of returned martian samples. Depending cm the 
size distribution of incoming particles and the role of natural 
resurfacing the quantity of “extramartian material could 
be significant. Analysis of samples will provide information 
on the meteoritic component as well as providing samples 
of unaltered added material. 

A. Banin discussed the mineral components of dust on 
Mars. Of most interest, he discussed recent laboratory 
analysis of martian analog materials in terms of their ability 
to reproduce the Viking lander labeled release experiment. 
Smectite clays satisfy all of the available data in terms of 
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surface composition and reactivity; interestingly, 
palagonites, which are a good spectral analog, do not fit 
the reactivity analysis and are therefore not a good analog. 

H. Moore discussed the global radar observations and 
implications for surface structure. Using a multicomponent 
system, with surface scattering analogous to lunar radar 
scattering, he was able to reproduce the observed radar 
variability by defining about two dozen surface units; each 
unit had a unique radar reflectivity and scattering ability, 
with the latter presumably related to some sort of surface 
or subsurface roughness. While the results may be 
nonunique, it is remarkable that the fit to the data was 
as good as it was. Further observations are required to 
better define the radar scattering from the martian surface. 

B. Jakosky described the variability of surface physical 
properties as seen in the remote-sensing data. Those parts 
of the surface that fit the general correlation of thermal 
inertia with radar reflectivity were identified, as were those 
parts that did not. Interestingly, the latter were located 
in discrete places rather than being scattered around the 
planet; these places were the canyon system and the crater 
Huygens (other locations probably exist, but the spatial 
extent of available observations was limited). These results 
have important implications for the ongoing surface 
evolution processes. 

J. Bell III summarized results of recent telescopic 
observations by himself and the group at University of 
Hawaii during the recent opposition of Mars. Observations 
were made at much higher spectral resolution than any 
previous ones and at spatial resolution comparable to the 
best seen from Earth. These new near-infrared reflectance 
data have a much higher signal-to-noise ratio than previous 
data and show spectral features not previously seen. Analysis 
of the data is ongoing, and interpretations of the results 
will appear shortly. 

DUST OBSERVATIONS FROM FUTURE 
SPACECRAFT MISSIONS 

The final session focused on future spacecraft 
observations: instrument status and capabilities, expected 
data, and possible collaborative activities related to dust 
studies. Three oral presentations were given. 

Zurek discussed the capabilities of the Pressure Modulator 
Infrared Radiometer (PMIRR), and Christensen elaborated 
on the expected observations to be made by the Thermal 
Emission Spectrometer (TES). Both experiments are now 
in construction and will be included in the Mars Observer 
payload. The capabilities of each instrument will be as good, 


or better than, those originally proposed. Of particular 
importance to dust studies and atmospheric modeling will 
be PMIRR’s determination of vertical profiles of 
temperature, dust content, and condensates as a function 
of atmospheric pressure. The TES will produce high- 
resolution compositional and thermal inertia maps of the 
martian surface, both of particular use for sediment 
transport and dust composition investigations. The planned 
Mars Observer mission will also allow spatial and seasonal 
variability of dust loading and transport to be determined. 

L. Esposito discussed the capabilities of the KRFM-ISM 
experiment on the Soviet PHOBOS mission. With spectral 
coverage out to 70 microns and observations covering a 
range of viewing and illumination angles, results from this 
instrument should further constrain models of dust size 
and distribution. 

All speakers lamented the deletion of the VIMS 
experiment from Mars Observer. Not only will important 
observations be lost, but intercomparisons with observations 
from the other Mars Observer and PHOBOS instruments 
will not be available. 

SUMMARY AND SUGGESTIONS 
FOR FUTURE WORK 

The benefits of collaborative research were stressed 
throughout the workshop. Inclusion of pertinent 
observations (such as dust properties, surface properties, 
patterns related to regional atmospheric circulation and 
sediment transport, timing and location of dust storms, 
and further detailed analysis of available Earth-based and 
spacecraft -derived data sets) with modeling efforts was 
thought to be particularly promising for future work. It 
was noted that the models of atmospheric transport and 
mixing are becoming increasingly sophisticated as they 
attempt to simulate more realistically the planetary 
boundary layer, dust and condensate aerosol physics, and 
the spatially three-dimensional nature of the atmosphere 
and of the atmosphere-surface interactions. Such models 
permit and require more detailed comparison with 
observations from which atmospheric winds and surface 
wind stress can be inferred. 

The following list enumerates some of the outstanding 
issues and how they may be addressed. 

1. The frequency of occurrence of the global and planet- 
encircling dust storms is still uncertain. Given the lack 
of NASA support for long-term monitoring of Mars from 
Earth (or near-Earth orbit), future observations are unlikely 
to refine adequately our current estimates. 
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2. The outstanding question regarding the evolution of 
a great dust storm is how it shuts off so quickly after onset. 
There are some interesting ideas, such as the depletion 
of mobile dust in the source region, >r negative feed-back 
mechanisms such as a fully nonlinear radiative-dynamic 
instability, but these need more work It would be desirable 
to find a single mechanism that can account for both the 
growth and decay phases of the life cycle of a great dust 
storm. 

3. With regard to the source(s) of airborne dust during 
a great dust storm, the question of whether the great storms 
can “exhaust” a region, only to have dust “resupplied” 
to that region within a very few years for the next storm, 
is an intriguing one. 

4 . We need better, more quantitative ideas of what 
happens at the air-ground interface to move dust. 
Eventually, we must have vertical probing of the planetary 
boundary layer (PBL) on Mars and a better description 
of its surface. Acquisition of appropriate terrestrial data 
could test significant aspects of the PBL models now in 
use. 

5. Is the zonally symmetric component of the martian 
atmospheric circulation largely an in viscid circulation? Or 
is it — at least at some latitudes — mechanically driven, and 
if so, through what mechanism (i.e., bi caking gravity waves)? 
Which circulation (inviscid or viscid) best represents the 
net transports of dust and volatiles? 

6. Finally, it is obvious that the N ASA/ Ames MGCM 
represents a major new tool for testing our understanding 
of the Mars atmosphere and of its interaction with the 
surface. Comparison of the MGCM results with both global 
and local data sets, like the distribution of bright and dark 
streaks, the Viking Lander meteorological data, and 
spacecraft observations of clouds and of temperatures are 
already underway. These types of studies will surely increase 
in number and in scope as people come to understand 
more fully the potential — and the limitations — of the 
MGCM. The model has shown sensitivity to specific 
dynamic-radiative feedbacks and to certain surface 
properties, like thermal inertia. The planned addition of 
a transport algorithm for dust and/or water vapor will 
improve significantly a vital means of resting our ideas about 
the cycles of dust, water, and carbon dioxide in the Mars 
atmosphere/surface system. 

The participants supported the idea of holding topical 
workshops similar to this one in the future. It was felt 
that the workshop format, allowing generous time for 
informal presentations and discussion, was very beneficial 
and maximized the opportunities for participants to 


exchange information and ideas. The feasibility of providing 
access to actual data sets (through the capabilities being 
developed by the Planetary Data System) in combination 
with presentations of analysis results (such as took place 
in this workshop) was discussed. Past experience with 
“presentation-oriented workshops” and “data analysis 
workshops” indicates the combination of the two formats 
would be conducive to even greater information exchange 
among the participants. During the next year, the possibility 
of organizing a future “Dust on Mars” workshop (probably 
in mid- to late 1990) will be explored. 
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THE MINERAL COMPONENTS OF DUST OH MARS; A. Banin, Dept, of 
Soil and Water Sciences, Hebrew University, 76100 Rehovot, Israel. 

The mineral components of the dust on Mars have not yet been 
directly analyzed. Up to the present, the major source of information 
on dust components is the various analyses, both in situ and by remote 
sensing, of the Martian soil, i.e. the fine, weathered, loosely packed 
aeolian material covering major areas of the planet. Evidence shows 
that this soil material has been thoroughly homogenized by the global 
dust storms characteristic of Mars, and it is plausible to assume that 
Mars dust composition is strongly correlated with the top-soil 
composition. This compositional similarity may be strengthened on 
Mars, where conditions over the last geological epochs did not favor 
rapid alteration and weathering of minerals either in the soil or in 
the dust. Various approaches have been used to constrain the 
mineralogical composition of Mars soil. These were based on chemical 
correspondence with the analyses of the soil done by the Viking 
Landers, simulation of the spectral-remote sensing observations of 
Mars, simulations of Viking Biology and other Viking Lander 
experiments, and various modelling efforts. On the basis of these 
studies, the following candidates have been suggested: 

1) Alumos ill cates : Smectite clays such as montmorillonite and the 
iron-rich nontroni te have been initially suggested as best satisfying 
spectral (1,2) and chemical (3,4) properties. Further strong support 
was obtained from simulations of the Viking Biology experimen ts ( 5-8 ) 
and the modelling of the weathering processes on Mars(9,10). Smectites 
do not however show the magnetic susceptibility measured on Mars. 
Palagonite, a poorly defined weathering product of basalts, is 
another candidate which satisfies spectral and compositional 
constraints ( 11-13). It fails, however, in the simulation of the 
reactivity shown in the Viking Biology experiments( 1 4, 1 5) • 

2) Iron Minerals:: Iron is a major chemical components in the Martian 
soil(3,4). Its minerals strongly affect the spectral characteristics 
of the soil in the visible range, lending the soil its characteristic 
reddish coloration. Reflectance spectra of the "bright" regions of 
Mars, believed to contain the more weathered components of the soil, 
show an intense though relatively featureless absorption edge from 
0.75um to the near ultraviolet ( 1 6 , 1 7 ) . Several mineral candidates 
have been shown to satisfy these spectral constraints including 
"amorphous" iron oxides ( 1 8 , 1 9 ) , palagonite ( 1 2 ) , corrundum doped with 
iron(20), jarosite(21 ) , cryptocrystalline hematite deposited in silica 
matrix( 22 ) and iron enriched smectites(7,8) . The major conclusion 
from the abundance of spectral analogs is that unambiguous 
identification of the iron mineral(s) in Mars soil, on the basis of 
reflectance only, is not possible at present. The spectral evidence 
does show, however, that iron in the weathered component of Mars soil 
is mostly present, in the oxidized form either in poorly crystallized 
clusters of oxy-hydroxy ferric iron, or as crystalline minerals, but 
in extremely small particle size range ( "nanocrystal3" ) . 

Smectite clsys chemically enriched with ferrous iron, which is 

then oxidized and polymerized in s_Ltu(23), 3how strong spectral 

analogy with Marti (Fig. 1). Compositionally , these preparations are 
full chemical analogs of the Mars soil with respect to iron oxide 
content. X-ray dif f ractometry , chemical extractions and electron 
micrographs show that iron in these preparations is in poorly 
crystalline or cryptocrystalline compound(s) and appears in particles 
of size range 1-10 nm(24). In these Mars soil analogs the 


PHECEOfM* p? l^y p* iVOT HIKED 


HfiE m INTENTtONAblY BLAH* 



16 


Minerals in Mars Dust 


Banin, A, 

crystallization of iron oxide minerals and the growth of their 
crystals is apparently prevented (a) by the presence of Si and A1 in 
the equilibrium solution, known to affect Fe oxide crystallization in 
soil environments(25) , (b) by the presence of an adsorbing surface of 
the smectite which immobilizes colloidal iron-oxyhydroxy particles and 
reduces iron concentration in solution, thus slowing down or 
preventing the precipitation and crystallization of pure iron 
minerals. Similar conditions can prevail in Mars soil, explaining its 
spectral properties and chemical reactivity. 

3) Carbonates: Although no spectral evidence for carbonates has been 
found yet on Mars (26) considerations of Mars climate evolution has 
led to suggestions that carbonate minerals should have been deposited 
on Mars(27; and shall, therefore, be present in the soil and dust on 
Mars. Recent simulations of the Viking Biology LR experiments using 
mixtures of the active ingredient (iron enriched smectite clays) with 
various carbonate minerals (Fig. 2) have shown that even low amounts 
(0. 5-1*0%) of added calcite reduce considerably the decomposition 
activity of the smectite and diminish the release of 14C0 2 , making the 

gher content of 
sible that the 
siderite rather 

than calcite, magnesite or dolomite. 

In summary, it is concluded that dust on Mars is likely to be a 
multicomponent mixture of weathered and non weathered minerals. 
Smectite clays are suggested as major components among the weathered 
minerals, assumed to be adsorbed and coated with amorphous iron 
oxyhydroxides and mixed with small amounts of better crystallized iron 
oxides ( f errihydri te , goethite and hematite) as separate phases, but 
having extremely small crystal sizes. As accessory minerals it is 
likely that the dust contains various sulfate minerals (Kieserite 
(MgSO-.HpO) and/or anhydrite (CaSO,)). If present, carbonates are 
likely to be at very low concentrations, although siderite (FeCO-,) may 
be present at Bomewhat higher concentrations. The chemical reactivity 
of this multicomponent system is complex to predict. It may be 
dominated by the ability of the smectite clays to adsorb volatiles due 
to their extensive surface area (up to 800 m2/g) and their catalytic 
reactivity in chemical and photochemical reactions. Furthermore, it 
is likely that the presence of the nanocrystalline iron oxy-hydroxy 
particles, with their large surface area, will increase the adsorption 
capacity of the system and enhance its reactivity. 
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Toon, O.B. et al. (1977) Icarus 30, 663 - 696 . 3 ) Clark, B.C. et al. 
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(1986) "Clays on Mars" in Clays and the Origin of Life (H. Hartman and 
A.G. Cairns-Smi th , eds.) pp lOo-l 1 5 , Cambridge Univ. Press. 9) 
Zolotov, M. et al. (1983) LPSC 14, 883-884. 10) Zolensky, M.E. et al. 
(1987) Mars Sample Return Workshop (Abst.) LPI/NASA/ JSC , 96-97. 11) 
Allen, C.C. et al. ( 1981 ) Icarus 45. 347-369. 12) Singer, R.B. (1982) 
JGR 87, 10,159-10,168. TJ7 Adams, J.B. et al. (1986) JGR 91, 
8098-8112. 14) Banin, A. and Margulies, L. (1983) Nature 305. 
523-526. 15) Banin, A. et al. (1988) LPSC 19, 27-28. 16) Singer, R.B. 
et al. (1979) JGR 84, 8415-8426. 17) McCord, T.B. et al. (1982) JGR 
87, 3021-3032. 18) Evans, D.L. and Adams, J.B. (1979) PLPSC 10, 
1829-1834. 19) Evans, D.L. and Adams, J.B. (1980) PLPSC 11, 757-763. 


mixture an unlikely candidate for Mars soil. A 
siderite (FeCO..) can be tolerated. Thus it 
carbonate component, if present in Mars soil, 


is 

may 


pos 

be 



17 


Minerals in Mars Dust 


Banir , A. 

20) Morris, R.V. el, al. (1983) LPSC 14, 526-527 21 ) I Burns R.G (1986) 

Nature 320, 55-56. 22) Morris, R.V. (1987) In, Mars — Sample _ Re tu rg^ 

Workshop LPI/NASA/ JSC Workshop Abstracts 96-97. 23) Qerstl » "* 

Banin, A. (1980) C lays Clay Min. 28, 335-345- 24) Ben 0 S . 1 ! 10 ” 0 ’ T T * " d 
A. Banin - unpublished data. 25) Schwertmann, U. (1988) in Iron in 
So ils and Clay Minerals (J.W. Stucki, B.A. Goodman and U. Schwertmann, 
eds. ) pp~ 267-308, D. Reidel, Dordrecht. 26) Singer, R.B. (1985) Mli 
Space Research 5(8) 59-68. 27) Kahn, R. (1985) Icar us 175-190. _ 

... . . ‘ • - NASA Exobiology and 


(This work was supported in part 
Systems Exploration Programs). 


by the 


Solar 



Fig. It Re f lec Uncp sjnctra of Mire and of 
ainecTlte ( montmorlllon 1 to ) cloy loaded with 
various amounts of iron in an acidic 
environment* Spectroscopically these Mara 
aoil analogs ehow reflectance similar to 
that of Mars soil at the wavelength range 
0.4-0.85um, due to oxy hydroxy-iron in 

poorly crystalline environments. 

Chemically, the hlphnst level of iron- 
addition cono is 1 1 tuie a a full chemical 
analog of Mars soil with respect to iron 
content. 



FI g . 2 1 The effect of carbonate minerals on 
The reactivity of a Mars soil analog in the 
almula tion of the Viking LR experiment on 
Mare. The anolog, smec ti te-mon tmorillon 1 te 
enriched to 10UX ion exchange capacity, 
containing in total 6,5X Fe^O^ shows 
reactivity similar to the Mara Soil (6,8), 
and lta decompoal tion capacity is taken as 
100X. Additions of carbonate decrease the 
decomposition as measured by 14C0- release. 
The affect of calcite is very significant 
even at low levels of additions, whereas 
Biderite can be present at much higher 
concentrations without interfering with the 
reaction , 


ORIGINAL PAGE IS 
OF POOR QUALITY 





BAROCLINIC-RADIATIVE INSTABILITY AND MARTIAN SUMMER 
POLAR CYCLONES; JR. Barnes, Department of Atmospheric Sciences, 
Oregon State University, Corvallis, OR 97331 

Some of the more interesting atmospheric phenomena observed by Viking 
were the polar cyclones (or so-called "spiral clouds"), seen at high latitudes in 
the summer northern hemisphere in several orbiter images (1 ,2). Based on the 
similarity in appearance to terrestrial extratropical cyclones, as well as the 
observation of sharp nearby temperature gradients by the Viking IRTM, these 
were identified as being manifestations of baroclinic instability in the s umm er 
polar atmosphere (1). The relatively small horizontal scales of the polar 
cyclones (~ 100-300 km) were noted as constituting a difference from terrestrial 
analogs. In contrast with the interpretation of baroclinic instability, the polar 
cyclones were alternatively proposed as being a result of radiative instability (2). 
This instability, as fundamentally investigated earlier (3), would be driven by 
atmospheric heating due to the presence of dust. Both the images and IRTM 
data indicated that a dust haze was present in the region of the polar cyclones 
(1 >2). A key difficulty with the radiative instability mechanism was that 
conditions of weak wind and weak baroclinicity (weak horizontal temperature 
gradients and vertical wind shears) were invoked as necessary to the 
development process (2). The IRTM thermal measurements indicated the 
presence of very strong baroclinicity in association with one of the storm 
systems (1). 

The apparently rather different dynamical explanations proposed for the 
summer polar cyclones have motivated some further studies, which are 
reported on in this paper. Specifically, combined baroclinic-radiative instability 
has been examined as a possible underlying dynamical mechanism for the 
polar cyclones, as well as a process of potentially broader relevance for the 
winter extratropical atmosphere. This instability has the basic character of 
baroclinic instability, as investigated previously for Mars (4), but is modified by 
the effects of dust heating (which by themselves lead to a pure radiative 
instability). 

A P-plane dynamical model like that used in linear baroclinic instability 
studies of the Martian atmosphere (4) has been employed for calculations of 
baroclinic-radiative instability. The effects of dust heating are simply 
parameterized in terms of the dust mixing ratio, which is determined 
according to a dust transport equation coupled to the dynamical model. The 
transport equation is the same (governing just the perturbation part of the dust 
field) as in earlier simulations of dust transport during a polar warming event 
(5). The coupled dynamical-dust transport model enables calculations of linear 
baroclinic-radiative instability to be made for relatively realistic basic states 
(zonal flows and dust distributions), with dissipation (radiative damping and 
surface friction). Without dust heating, the model yields unmodified baroclinic 
instability as examined previously (4). With dust heating, but without 
baroclinicity (that is, without vertical shear of the basic zonal flow), the model 
yields pure radiative instability fundamentally like that investigated earlier 
(2,3). In the presence of both dust heating and baroclinicity, combined 
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baroclinic-radiative instability can occur. It can be anticipated, on the basis of 
the relatively low growth rates for pure radiative instability (2), that the basic 
character of the combined instability will be that of barochmc instability (for 
reasonably large barodinicity and dust heating and dust gradients that are not 
overly strong). It will be modified, though, by the effects of the dust heating. 
This modification may be especially important at the relatively short scales of 
the summer polar cyclones. 

Results from preliminary calculations with the p-plane dynamical- dust 
model have revealed several aspects of radiative instability that were not 
apparent from the earlier studies. The first is that the growth rate spectrum 
for purely radiative instability is flat, the instability having virtually no scale 
dependence. The existence of a nearly flat growth rate spectrum in the 
presence of a significant P-effect was actually indicated by the previous simple 
model studies (3). To further investigate the scale dependence of radiative 
instability additional calculations were performed with another simplified 
model, a two-level dynamical-dust model. In this model the radiative growth 
rate spectrum falls off' at the longest scales, with the fall-off becoming more 
pronounced as the strength of the instability (proportional to the dust heating 

and the dust gradient ) is increased. The spectrum becomes flatter as the P* 
effect increases. The growth spectrum in the more realistic (multi-layer) 
model, however, is essentially completely flat both with and without a P-effect. 
The key differences between the simplified models and the present more 
realistic one appear to be associated with vertical structure. 

The other important aspect of the radiative instability problem that is clearly 
seen in the multi-layer model results is that the vertical structure of the 
unstable radiative modes is tied very closely to the vertical profile of the basic 
state dust distribution. The modes have greatest amplitudes where the (basic 
state) vertical dust gradient is largest. In contrast, baroclinically unstable 
modes have their largest amplitudes at relatively low levels in essentially all 
cases. A consequence of this is that a shallow dust distribution (having a 
maximum vertical gradient in the lowest scale height) is required in order for 
the effects of dust heating to be significant for the baroclinic modes. 
Additionally, the dust heating (and vertical dust gradient) must be fairly strong 
and/or the barodinicity fairly weak in order for the effects to be important. If 
these conditions prevail then dust heating and radiative instability can play a 
significant role in the dynamics of atmospheric disturbances. Northern 
summer middle and high latitudes on Mars should indeed be characterized by 
relatively strong dust heating and weak baroclinirity, as well as relatively 
shallow dustiness, so that dust heating and radiative instability may be of at 
least some importance for the dynamics of summer polar cyclones. 
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PRELIMINARY RESULTS AND OVERVIEW OF SPECTRAL REMOTE 
SENSING OBSERVATIONS OF MARS DURING THE 1988 OPPOSITION; James 
F. Bell III, Thomas B. McCord, Paul G. Lucey, and Diana L. Blaney (Planetary 
Geosciences Division, HIG, Univ. Hawaii, Honolulu, 96822) 

Numerous research and observational studies of Mars are being carried out at 
the University of Hawaii Planetary Geosciences Division (PGD) during the excellent 
1988 opposition. A combination of high angular diameter, new high spectral resolution 
telescopic instrumentation, and Mauna Kea Observatory’s low atmospheric water 
content provide a unique opportunity to resolve some of the presently unknown or 
poorly understood aspects of martian surface chemistry and mineralogy. 

Studies underway include (1) Moderate spectral resolution (AA/A~1.5%) point 
spectroscopy of 200-300 km size regions in the 0.3-2.5 /im wavelength range. Goals 
include characterization of specific Fe 2+ and Fe 3+ mineralogies in the 0. 3-1.0 /xm region, 
analysis of H 2 0 and C0 2 frost distributions and/or concentrations in the 1.0-2. 1 M m 
region, and identification of specific cations associated with previously observed yet 
uncharacterized clay absorptions in the 2.2-2.5 Mm region [1,2]. (2) High spectral 
resolution (AA/A~0.3% j imaging spectroscopy of (a) 150-200 km size regions at 0.3-1.0 
/xm to identify suites of Ft 2+ and Ft ^ absorptions and thus place severe constraints 
(based on laboratory studies [3,4,5]) on the composition and distribution of crystalline 
Fe mineralogies in the Mars dust and/or soil, (b) 300-500 km size regions at 2. 2-2. 5 pm 
to identify suites of clay mineral absorptions and the specific cations responsible for the 
cation — OH absorption features seen at these wavelengths, (c) 300-500 km size regions 
at 3.0-5.2 M m to search for possible absorption suites characteristic of carbonates, 
nitrates, or sulfates postulated to exist on the surface of Mars and shown to be 
observable from ground based observations if existing in abundances >~3-5% [6]. These 
imaging data will be incorporated into a set of equal-area projection maps depicting the 
distribution of the various inferred mineralogies as a function of martian time of day as 
well as L,. (3) Discrete filter imaging in the 1. 2-2.6 fim region using an InSb linescan 
array [7] to produce maps similar in form to those discussed above of the distribution 
and variations in clay mineralogies at 300-500 km resolution and H 2 0 and C0 2 frost 
deposits at sub-hemispheric scales. (4) Moderate spectral resolution (AA/A~1.6%) 
infrared spectroscopy in the 7-11 fim region to characterize the position of the 8/im 
martian Christ iansen emission feature. 

These studies will be discussed in greater detail during this workshop 
presentation and preliminary results dealing with moderate resolution point spectra in 
the 0.3-1.0 ptm region will be presented. 

References: (1) McCord, T.3. ct of. (1982) J. Geophys. Res., 87, 3021-3032. (2) Singer, et a L (1979) J. 
Geophys. Res., 84, 8415-8428. (J) Hunt, G.R. and J.W. Salisbury (1970) Mod. Geol., 1, 283-300. (4) Ibid 
(1971) Mod. Geol., £, 195-205. (5) Sherman, D.M. et al. (1982) J. Geophys. Res., 87, 10169-10180. (6) 
Walsh, P.A. and T.B. McCord (1987) EOS, Trans. Amer. Geophys. U., 68, 1341. (7) Lucey, P.G. et al. 
(1987) B.A.A.S., 19, 843. 
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MARS OBSERVER VIMS AND PHOBOS KRFM-ISM OBSERVATIONS OF MARS DUST; 
L.W. Esposito, LASP/Univ. of Colorado 

The visible and near-IR is a fruitful region to observe dust because of 
its substantial opacity at these wavelengths, compositional implications from 
near-IR bands, and the possibilities of volatiles coating the dust particles 
Because of the varying strengths of C0 2 bands in this region, the total abun- 
dance of dust and its vertical distribution can be inferred from comparing 
the band strengths and shapes. This method is useful even at nadir viewing 
and thus it can measure larger columns of dust (such as during dust storms) 
that cannot be probed by limb viewing. The Mars Observer orbit and mission 
length allow global and seasonal variation in dust to be determined. The 
high spatial resolution of VIMS can identify heterogeneity and allow small 
storms to be studied. The KRFM-ISM experiment on PHOBOS has lower spatial 
and spectral resolution, but has the advantage of observing at a range of 
viewing and sun angles (further, it is now en route to Mars). KRFM bands up 
to 70^ will further constrain dust size and distribution. If VIMS or a 
scaled-down spectrometer is a part of the Mars Observer science payload, 
valuable intercomparisons are expected. 
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AN ASSESSMENT OF THE METEORIC AND METEORITIC CONTRIBUTIONS 
TO DUST ON MARS; G.J. Flynn 1 , and D.S. McKay 2 , (1) Dept, of Physics, SUNY-Plattsburgh, 
Plattsburgh, NY 12901, (2) SN14, NASA Johnson Space Center, Houston, TX 77058 

Four types of meteoric and/or meteoritic material should be found on Mars: 1) micrometeorites, 
many of which will survive atmospheric deceleration unmelted, which should fall relatively uniformly over 
the planet’s surface, 2) ablation products from larger meteors and meteorites which ablate, break up 
and/or burn up in the Mars atmosphere, 3) debris from large, crater forming objects, which, by analogy to 
terrestrial and lunar impact events, will be concentrated in the crater ejecta blankets (except for rare, large 
events, such as the proposed C-T event on earth, which can distribute debris on a planetary scale), and 
4) debris from the early, intense bombardment, which, in many areas of the planet, may now be 
incorporated into rocks by geologic processes subsequent to the intense bombardment era. 

To estimate the extent of meteoritic addition to indigenous Martian material, the meteoritic flux on 
Mars must be known. Hartmann et al. (1981) estimate that the overall flux of large, crater producing 
objects at Mars is twice that for the Moon and 1 .33 that for Earth. For small particles (< 1 00 micrometers 
in diameter), whose orbital evolution is dominated by Poynting-Robertson (PR) drag (Dohnanyi, 1978), we 
estimate the flux at Mars from the measured flux at Earth. The smaller Martian gravitational enhancement 
as well as the decrease in the spatial density of interplanetary dust with increasing heliocentric distance 
(inferred from the fall off of zodiacal light intensity) should reduce the flux of small particles at Mars to 
about 0.33 times the flux at Earth. Because of the smaller planetary cross-section the total infalling mass 
at Mars is then estimated to be 0.09 times the infalling mass in the micrometeorite size range at Earth. 

INFALL ASSESSMENT: At Earth the current annual infall of meteoric and meteoritic material hitting 
the top of the atmosphere is estimated, from satellite, radio and photographic meteor measurements, at 
16,000 tons, with 12,400 tons of that concentrated in the particles in the 10*° to 10* 2 gram mass range 
(Hughes, 1978). Ir accumulation rates in sea sediments suggest a relatively constant meteoritic and 
meteoric infall on earth over the past 67 million years (Kyte, 1 986), though Kyte infers a long-term infall 
rate of 77,000 + 25,000 tons per year, about five times the current infall estimated by Hughes (1978). The 
10' 6 to 10* 2 gram particles (from 100 to 3000 micrometers in diameter for density 1 gm/cm J ), which 
contribute 75% of the extraterrestrial material, are somewhat larger than the size whose orbital evolution 
is dominated by PR drag. The Mars flux of these particles is thus likely to be between that which we 
calculate for PR drag dominated particles and that suggested by Hartmann et. al. (1 981 ) for the larger 
objects. Taking the 0.09 reduction factor we calculate for the smaller particles as a lower limit for the 
Mars/Earth infall ratio for 10*° to 10* 2 gram particles, and using the lower value of the Earth infall rate from 
Hughes (1 978), a minimum of approximately 1 1 00 tons of meteoritic material in the 1 0* b to 1 0* 2 gram 
mass range would be deposited on the top of the Martian atmosphere annually. If this mass reached the 
surface as solid material, it would be sufficient to cover the planet with a 3 to 4 cm thick layer over the age 
of the solar system. This micrometeoritic infall amounts to about 0.7 cm/10 years. Substantially larger 
values are possible if either the Hartmann et al. Mars/Earth flux ratio or the higher Kyte earth flux is 
adopted, with an accumulation of about 14 cm/10 9 years resulting if both values are used. 

DIRECT COLLECTION: On earth, particles in the 10* 6 to 10* 2 gram mass range are normally melted or 
volatilized on atmospheric entry, with only particles < 100 micrometers in diameter surviving entry 
relatively unaltered (Brownlee, 1985). However Mars, because of its low gravitational acceleration 
combined with sufficient atmospheric density to provide deceleration, is probably the most favorable site in 
the solar system for unaltered survival (and thus collection) ol micrometeorites. We have made new 
calculations of the interactions of these micrometeorites with the Martian atmosphere using the 
micrometeorite deceleration model first used by Whipple (1950) and the upper atmospheric density 
profile for Mars derived from Viking entry measurements (Seiff and Kirk, 1976). Particles of the same 
density, shape, and thermal characteristics entering at velocities near the Martian escape velocity are 
heated to a peak temperature only half that experienced on Earth entry at earth escape velocity. 

Although most particles larger than 75 to 100 micrometers in diameter are melted on Earth atmosphere 
entry, the cutoff size for melting on Mars entry would be about 13 times larger, or 1000-1300 micrometers 
in diameter. 

The concentration of micrometeorites in Martian soils depends on the fraction of them which survive 
atmospheric entry as solid material (either unmelted, as are the micrometeorites collected from the 
Earth’s stratosphere, or melted, as are the chondritic spheres collected from the sea sediments - see 
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review by Brownlee (1985)), as well as the Martian regolith depth, and the total accumulation time. For 
illustrative purposes, we assume that all of the micrometeorites smaller than 3000 micrometers survive 
atmospheric entry in some form as solid material, the Martian regolith has an average depth of 10 m 
(consistent with Arvidson’s (1986) estimate of a planetary erosion rate of meters/10’ years, based on 
crater preservation), and that micrometeorites have been accumulating at the present estimated rate for 
the past 4 x 1 0 a years. For these conditions, a 1 0 gram average soil sample would have been mixed with 
about 5000 micrometeorites greater than 100 micrometers in diameter and 10 micrometeorites greater 
than 800 micrometers in diameter. Micrometeorites in this size range which have survived atmospheric 
entry are normally not collected at earth. Since these larger particles from Mars orbital distance are likely 
to sample different sources than the smaller micrometeorites collected at 1 a.u. in the cosmic dust 
sampling program on earth (Flynn, 1988; Zook and McKay 1986), returned Mars soils may provide a 
unique resource for micrometeorite analysis. 

MICROMETEORITES AS A TOOL: The degree to which Martian weathering will break down 
micrometeorites is unknown. However, since unmelted 1000 micrometer fragments of meteoritic material 
have been identified in the 2.3 million year old Pliocene layer on Earth (Brownlee et al., 1982), 
micrometeorites added to Martian regolith after atmospheric deceleration below hypervelocities may still 
be identifiable by petrographic or chemical means. Assuming that micrometeorites could somehow be 
identified in returned soil samples, this addition of micrometeorite material to the uppermost Martian 
regolith at a constant rate could conceivably provide a powerful tool for tracking rates of erosion, 
deposition, and weathering. On Mars sample return missions, an attempt should be made to collect soils 
from different geologic sites (catch basins, lag surfaces, flat high plains, valley bottoms, etc.) so as to 
provide a variety of soils of different sedimentary environments. One of the important differences among 
these environments might be the proportion of petrographically or chemically identifiable 
micrometeorites mixed into the soil. 


COLLECTING SITES FOR MICROMETEORITES: Martian surface processes (weathering and wind 
erosion, transport, and deposition) may fractionate the dust by size, density or composition providing 
regions of increased local concentration, suggesting even more suitable sites for micrometeorite 
sampling than the average soil. These sites may include placer catch basins or lag surfaces which may 
accumulate high density micrometeorites or their derived and altered minerals. Conversely, low density 
micrometeorites may be wind segregated along with finer Martian dust and may constitute a relatively 
coarse-grained component of that dust at its deposition sites. By analogy with Antarctica, meteorites of all 
size ranges may be relatively concentrated in Martian polar regions, although the concentration 
mechanisms may be different. 

CHEMICAL SIGNATURES: The possibility that detectable micrometeorites and their remains can be 
found in the Martian soils depends on the relative rates of infall, weathering and alteration, transportation, 
and mixing. These rates are not yet known reliably enough to allow us to predict with certainty whether 
identifiable micrometeorites will be found. While they may be relatively quickly destroyed by Martian 
weathering, the chemical signature, particularly siderophiles and volatiles, may persist in the soils, as 
they have in the lunar regolith (Anders, et al., 1973). If the regolith production rate on Mars is of order 1 
meter/10 9 years, our calculated infall rate would give rise to a mean concentration of about 1% 
micrometeoritic material in the regolith. 

In the lunar case the composition of the non-indigenous material was taken as the residual after 
subtracting rock composition from soil composition. Two meteoritic components were detected. In 
mature soils the residual has a trace element composition consistent with the addition of 1 .5% Cl 
meteoritic material, attributed to the long term micrometeorite infall (Anders et al., 1973). Less mature 
highland soils also show a second component, characterized by fractionated siderophile content and low 
volatiles, attributed to ancient bombardment (Anders et al., 1973). 

Boslough (1988) and Clark and Baird (1979), applying the subtraction method to the Viking chemical 
data, suggest the Mars regolith can be fit by a mixture of 40% Cl meteorite and 60% planetary rock 
fragments. Boslough (1988) suggests the meteoritic component is ancient. But it could equally well be 
the micrometeorite component, which dominates the ancient component in lunar soils. If either the 
planetary regolith production rate is significantly lower than the 1 meter/10 9 year value we have used for 
illustrative purposes or if the micrometeorite infall is higher than we calculate using conservative values 
for the earth flux and the infall ratio, the fraction of meteoritic material would increase from the 1% we 
calculated. Increasing the earth infall rate to the value reported by Kyte (1986) and the Mars/Earth infall 
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ratio to the Hartmann et al. (1981) value, and reducing the Mars erosion rate to one-third of the 1 
meter/10® years would produce a micrometeorite concentration in the Martian regolith of about 40%. 

Thus a Martian regolith containing a substantial fraction of meteoritic material cannot be ruled out. On 
earth the Ni/Fe ratio and the Ir abundance have proven to be diagnostic indicators of meteoritic material, 
since Ir and Ni are enriched in Cl meteorites but depleted in crustal materials. Direct measurement of the 
Ni and/or Ir abundances in the Mars regolith should help determine the meteoritic content of the regolith. 

The micrometeorites collected at Earth, because of their small size, contain large concentrationss of 
solar wind implanted ions, particularly the noble gases (Hudson et al., 1981; Rajan et al., 1977) and likely 
hydrogen. These solar wind species, carried into the Martian regolith by the micrometeorites, may also 
perturb the bulk chemical abundances. 

MARTIAN AGGLUTINATES: If, as we calculate, micrometeorites are all slowed down by the Martian 
atmosphere, and assuming that most lunar agglutinates are made by micrometeorite impacts, no 
analogous Martian agglutinates would be expected (unless there were an era in which the atmosphere 
was considerably less dense than at present). However, many types of impact glasses would be expected 
from larger impacts, and some of these glasses may resemble lunar agglutinates in some respects. 

MARTIAN SOIL MATURITY: Gault and Baldwin (1970) have estimated a minimum impact crater 
size of 50 meters, taking into account fragmentation and ablation of the incoming projectiles as well as 
atmospheric deceleration. The smallest craters noted in Viking orbiter images are about 100 meters in 
diameter (Blasius, 1976), but smaller craters beyond the resolution limit of the photographs may still be 
present. Dycus (1 969) predicts that projectiles as small as 10 gm would still form craters. However, 
craters too small to be seen from the orbiter are not apparent in Viking lander images. Impact gardening 
associated with the 50 meter and larger craters predicted by Gault and Baldwin (1970) would determine 
regolith turnover rates and cause comminution of rocks into soils. The addition of micrometeorites would 
affect the petrology and chemistry of Martian soil. Weathering and sedimentary processes on Mars would 
also process the regolith components. The overall effect would be to make an exceedingly complex 
regolith. A new maturation scale will be necessary for Martian regolith. This scale will have to include 
terms which reflect (1) impact reworking, (2) addition of micrometeorites, and (3) Martian surface 
weathering and alteration. F or example, if concentration mechanisms can be factored out, the abundance 
of micrometeorites (identified petrographically or chemically) in a soil layer might be directly related to its 
near-surface exposure time in a manner analogous to the abundance of agglutinates in lunar soils. In 
addition to soil evolution through maturation, physical mixing of soils of differing maturities should be 
common. 

CONCLUSIONS: The addition of meteoritic material to the Martian regolith could significantly perturb 
the chemical abundances in the soils, particularly the abundances of volatile and siderophile elements 
which are abundant in Cl meteorites but depleted in crustal materials, and of noble gases (and possibly 
hydrogen) which are implanted in the micrometeorites during space exposure and carried into the soils 
with the particles. The first returned soil samples from Mars should provide the opportunity for recovery and 
analysis of unaltered microrneteorites larger than any sampled on earth, assessment of the magnitude of 
the meteoritic component, end possibly an estimate of the rate of erosion and regolith production on the 
planet. This micrometeorite population may be quite different from the population sampled at 1 a.u. The 
extent of regolith gardening, small crater production and agglutinate production, if any, will also provide 
clues to the evolution of the Martian atmospheric density over time. 
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fhotographs of the Martian surface reveal the coexistence of a red 
powdery sediment and larger pebble and octoble-like elements. For the 
atmospheric density of mars, it is unlikely that these large elements are 
moveable by the wind. Nonerodible elements have a significant affect on dust 
production on the earth and presumably on Mars. 

The effect of nonerodible roughness elements on a surface of erodible 
particles is to shelter part of that surface by absorption of part of the wind 
momentum flux. In cur wind tunnel experiments dene in Boulder, 00, spheres 
too large to erode were half buried and appeared as hemispheres above a flat 
surface formed of erodible spheres. The effect of the nonerodible spheres on 
the saltation flux of erodible spheres was to increase the apparent threshold 
friction velocity for wind erosion. Above the threshold velocity, the 
saltation flux behaved as an Owen function as follows: G = A U* (U * 2 - U* t 2 ) 

where G is mass flux, A is a constant, U* is friction velocity, and U* t is 
apparent threshold friction velocity. 

Comparison of our results with those of the momentum stress partitioning 
experiments of Marshall (1971) showed that for particles having a small 
threshold velocity, the nonerodible roughness elements afforded less 
protection than predicted by Marshall's momentum partitioning relationships. 

We explain this as reflecting the inertia of the measuring devices: the fine 

sand used in our experiments furnished hi^ily sensitive measurements of 
momentum partitioning; the stress plates used in the Marshall experiment were 
apparently less sensitive measuring systems. 

The apparent threshold friction velocity may be obtained by a knowledge 
of size distributions of aggregates and the physical state of the soil. Fran 
this, expected wind erosion may be calculated from a knowledge of the wind 
probability density function. 


POSSIBLE INTENSE VORTICES AND THE POTENTIAL FOR DUST AND 
SAND TRANSPORT ON MARS; J.A. Grant and IML Schultz, Geological Sciences, 
Brown University, Providence, RI 02912 

INTRODUCTION: Characteristics of distinct, dark, ephemeral filamentary linea- 
tions in certain locations on Mars and their occurrence only during specific seasonal 
atmospheric conditions indicate they result from intense atmospheric vortices (1). 
Similar markings have been interpreted as linear seif dunes (2) or joint patterns (3), but 
these explanations are inconsistent with all observations. If atmospheric vortices of up 
to tornadic intensity do exist on Mars they should be capable of lofting and transporting 
significant amounts of dust and sand-sized material (4). 

DESCRIPTION: Individual lineations have total lengths ranging from 2 km to over 
75 km and widths up lo a few hundreds of meters (Fig. 1). They have well defined, 
straight to curvilinear plans that commonly cross each other and traverse considerable 
relief. At the highest resolution available (-75 m/pixel) no relief on the lineations has 
been identified. While the majority are found within several midlatitude areas in the 
southern hemisphere, their occurrence is quite widespread (but rare near the equator). 
The lineations develop with a predominantly E-W or NE-SW orientation on stripped 
intercrater plains and the floors of some large craters during the late summer: they then 
are rapidly modified and become undetectable by midfall (1). Lineations appear to 
redevelop in the same broad areas from year to year, but the location of individual linea- 
tions within these areas changes. 

DISCUSSION: As summarized in (1), the lineations are not thought to be the ex- 
pressions of joints or dunes (2,3) because the markings: (A) are highly variable in occur- 
rence over short time scales; (B) traverse landforms, each other, and local structural 
trends, without deflect on or interuption (Fig. 1); (C) change position from year to year; 
(D) are insensitive to broad scale structural trends; (E) are orders of magnitude larger 
than most terrestrial joints; (F) have grossly similar orientations over broad areas; and 
(G) cross materials whose competence should be too low for joint expression (i.e. ejecta 
facies). The occurrence and characteristics of all of the lineations indicate that they 
result from surface scour and re-deposition of material during the passage of intense 
atmospheric vortices (i). Considerable fine-grained material would be lofted and carried 
along by these vortices thereby leaving a stripped, coarser grained surface lag. Coarser 
sand-sized material also should be lifted (4), but would experience more limited 
transport. The resultant bands of relatively coarse-grained material remaining behind the 
intense vortices would be dark compared to the surrounding, unaffected areas (5). 

Although the occu rence of dust devils on Mars has been confirmed (6) it is unlikely 
that they are responsible for the formation of the lineations. First, the size (length and 
width) and the some of the characteristics of the martian lineations are inconsistent with 
those of martian and terrestrial dust devils (1,6,7). Second, the observed martian dust 
devils did not produce surface lineations. And third, atmospheric conditions present 
during formation of the lineations are not optimal for dust devil formation. 

Martian lineations begin to form in late summer following dust storm activity (1) 
when the clearing atmosphere allows development of a deep convective layer (8). The 
occurrence of this deep convective layer in combination with unstable atmospheric con- 
ditions (1) and the omet of seasonal baroclinic wave passage (9-11) establish conditions 
similar to those on the Earth favoring tornadogenesis (12). Consequently, strong at- 
mospheric vortices of possible tornadic-intensity might also form on Mars (1). While 
terrestrial tornadoes aio driven by latent heat released by water condensation, martian 
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vortices would be driven primarily by extended dry convective uplift of unstable surface 
air in the deep convective layer. Only a minor contribution at most would be made by 
latent heat release due to the paucity of water in the martian atmosphere (13,14). 
Because these intense vortices are controlled by midlevel rather than near-surface at- 
mospheric processes (12), they would be less sensitive to topography than dust devils and 
therefore more likely to produce groundtracks which traverse landforms and have gaps 
unrelated to topography. Finally, terrestrial tornadic vortices generate surface markings 
that resemble the martian lineations (12). 

The apparent absence of lineations on Mars during the spring when similar at- 
mospheric conditions exist may be related to the scarcity of images in those locations 
where they are expected to form (1). Complete destruction of the lineations by midfall 
is likely the result of burial by dust from atmospheric fallout including processes as- 
sociated with the expanding seasonal polar cap (1). 

IMPLICATIONS: From the distribution of lineations we estimate that they form 
seasonally in areas covering roughly 1.7 X 10 s km 2 , or about 1.2% of the surface of Mars 
with the largest fraction (>82%) occcurring between 30°S and 65°S. The average linea- 
tion density (55 to 60 per thousand square kilometers) in the latter areas combined with 
their average width (-200 m) and length (-10 km) represents about 2 X 10 s km 2 total 
area and permits an estimate of the seasonal eolian transport by this processes. A value 
for the average erosion resulting from the passage of a martian intense vortex of 1 mm 
was utilized in these calculations since measurements of the load carried by intense vor- 
tices on the Earth are not readily available. The area of each track multiplied by the 
amount of erosion gives a total average dust and sediment load of 2 X 10 3 m 3 entrained 
by each vortex over a distance of 10 km with the maximum mass lofted about 4.0 X 10 3 
kg/vortex (for material density of -2.2 g/cm 3 ). For comparison, the load carried by 
martian dust devils has been estimated to be 3 X 10 3 kg each based on their optical depth 
(6). The value given for the load carried by each martian intense vortex is the sum for 
the entire path-length travelled. The load carried by a single intense vortex at any mo- 
ment should be less than the maximum mass predicted because much of the entrained 
fine grained material (<5-10 microns) will be lost through the top of the vortex while 
some of the coarser material will simply be redeposited along the track. For the total 
area seasonally covered by lineations as much as 0.18 km 3 of dust and sediment could be 
lofted over a short period during late summer. 

The fine grained fraction lofted by these vortices may serve to prolong the decay of 
large martian dust storms or possibly contribute to their redevelopment. If intense vor- 
tices also develop during the southern hemisphere spring as analogues to seasonal tor- 
nadogenesis on the Earth, then they could play an important role in the initial develop- 
ment of certain large dust storms. Dust lofted in the Hellespontus region could contribute 
directly to dust storm development while dust raised in midlatitudes further to the west 
may have to drift north towards Solis planum before contributing to development. Lofted 
materials of all sizes would be transported away from source areas where the intense 
vortices develop. While fine grained material could contribute directly to the formation 
of polar layered terrains, coarser material could accumulate on crater floors and in the 
circumpolar dune fields, subsequently reworked by other winds to form dunes. 

Although the seasonal erosion due to intense vortices on Mars will be small, the 
cumulative yearly erosion could be significant. Such vortices may be responsible for 
much of the stripped appearance of the northern martian plains, but the lower albedo 
produced by removal of the finer material may preclude detection of the lineations there. 
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Fig. 1. Filamenta y lineations located at 91°W, 48.5°S just over 1000 km west of 
Argyre basin. Lineations trend E-W across the image and cross both crater walls and 
large N-S oriented frac tures. Viking image 523A34 
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COMPARISON OF MARTIAN AEOLIAN FEATURES AND 
RESULTS FROM THE GLOBAL CIRCULATION MODEL; R. Greeley, 

rrtvtiL uFo a °{ Geo l° gy ’ Anzona State University, Tempe, AZ 85287-1404, and 
J.B. Pollack, NASA-Ames Research Center, Moffett Field, CA 94035. 

Aeolian (wind) processes link the atmosphere of a planet with its surface. Wind is 
apparently the dominant agent of surface modification on Mars today, as evidenced by frequent 
dust storms. Features attributed to aeolian processes include dunes, yardangs, grooves, deflation 
pits, and albedo patterns that are time-variable. In addition, vast regions appear to be mantled 
with sediments presumably deposited from the atmosphere. Because of the common occurrence 
°f »u° i- an features on Mars and the frequent occurrence of dust storms, understanding the nature 
ot the features and of the interaction between the atmosphere and lithosphere are critical for 
several aspects of martian science: (1) the presence or absence of active and inactive dust mantles 
may affect the interpretation of some remote sensing data for Mars, (2) the location and 
occuirence of aeolian deposits will strongly influence the selection of potential landing and 
sample return sites, and (3) interpretation of martian surface evolution must include the 
identification and separation of material units formed by a wide variety of processes (fluvial 
penglacial, volcanic, etc.). ’ 

w S ^ a o C c e / e J atUre ! C j n ** used to “fier wind directions and strengths. For example, 

Ward et al. (1985) described yardangs oriented with the prevailing winds at the time of 
onnation. However, they noted several "sets" of yardang orientations and suggested that there 
had been shifts in wind patterns with time. Wind streaks probably constitute the best "wind 
vanes . As described by Thomas (1982) and others, many wind streaks change with time and 
may reflect seasonal shifts in wind patterns. Wind streaks and the interpretation of wind patterns 
^ JSf relatively well documented by Sagan et al. (1973), Thomas et al. (1981), Veverka et 
al. (1977), Greeley et al. (1977), and others. Incorporation of Viking IRTM data have enabled 
assessment of sediment sources, transport paths, and deposition sites of windblown sediments 
(Christensen, 1986; Lee, 1986). Dunes occur in high to mid latitudes (Breed et al., 1979- 
fhomas, 1982) and in small fields in equatorial regions. Although most are crescentic and 
barchan types (Breed et al., 1979; Tsoar et al., 1979) indicating a unidirectional wind regime 
reversing winds were also suggested by Tsoar et al. (1979). Cutts and Smith (1973) and 
Lancaster and Greeley (1987) identified reversing and star-like dunes in some areas, suggesting 
seasonal changes in wind direction. However, the relationship between dune patterns and 
formative winds remains uncertain and there is poor agreement between the alignment patterns of 
die dunes and some models of wind patterns (Ward and Doyle, 1983). This may suggest that the 
dunes formed when martian wind regimes were different. Breed et al. (1979) suggest that some 
North Polar dunes are eroding today, or being modified under long-term changes in wind 
regimesjTsoar et al (1979) and Ward and Doyle (1983) suggest that some dunes £e currently 
active. Thomas (1981) suggests southern hemisphere dunefields are aligned with current winds. 

Important questions remain about the nature of dunes and other aeolian features on Mars 
and their relationship to past and present winds. The Global Circulation Model (GCM) 
winom a Pollack , et ^ 1S a complex simulation of the martian atmosphere based on Toon et 
at. (198U). Among the parameter predictions made by this simulation are predictions of wind 
surface shear stress. A program running on the Cray 2 computer allows predictions to be made 
for the wind surface shear stress as a function of atmospheric temperature and pressure, location 
on Mars, and other parameters such as topography, atmospheric dust-loading, and local surface 
roughness. Output from the program includes maps with wind vectors and strengths for stated 
conditions of martian season and climate. The magnitude, position, and azimuth of the GCM 
stress vectors are read from tape and are projected to a mercator base map. Wind streak and other 
aeolian data are digitized from published maps. The position and azimuth of these data are then 
calculated and projected to match the base map. Map symbols representative of aeolian features 
(bright streak, dark streak, yardang, etc.) are then drawn. 
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Figure 1 shows the GCM surface stress data and digitized bright streak data plotted cm a 
mercator base map. The GCM data are arranged in 7.2° latitude by 9.0 longitude bins. The 
wind barbs indicate the direction and magnitude of the surface stress, and represent averages ot 
discrete stress values produced for the simulated period: Ls = 281.02 - 286.73, sampled at . 
hour increments. The bright arrows indicate downwind direction as revealed by bright streak 
orientation. The aeolian data used for Figure 1 are bright streak azimuths digitized from a globa 
compilation of aeolian features by Ward et al., 1985. The results from these comparisons will 
enable the following questions to be addressed: (1) Is there an alignment of aeolian features with 
predicted wind directions?, (2) Do erosional features occur in areas of P^cted sttong [winds/, 
(3) Do aeolian depositional features occur in areas of low winds?, and (4) Do theGCM transport 
paths correspond to trends suggested by sediment transportation derived from IRTM studies such 
as Christensen (1983, 1986) and Lee (1986). 
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Figure 1. GLOBAL CIRCULATION MODEL AND BRIGHT STREAK DATA FOR THE THARSIS REGION 


GCM Parameters: 

Ls (deg) 281.02 - 286.73 
Sampling increment - 1.5 hours 
Stress averages for the simulation 
are displayed. 


Wind barb explanation: 

stress units * N/m 2 x 10‘ 4 
short dash - 5 units 
long dash - 10 units 

triangle - 50 units 


Bright streak symbols indicate 
inferred down-wind direction, 
data from Ward et al. t 1985. 

(white) N 



MARTIAN AEOLIAN FEATURES 
Greeley, R. et al. 


33 


SIMULATIONS OF THE MARTIAN BOUNDARY LAYER: FACTORS CONTROLLING THE 
BEHAVIOR OF THE SURFACE STRESS; R.M. Haberle, J.B. Pollack, and J. 

Schaeffer, Space Science Division, NASA/Ames Research Center, Moffett Field, 
CA, 94035. 

We have begun to explore the properties of the Martian atmospheric boundary 
layer using a one-dimensional (height) time-dependent second-order 
turbulence closure motel that accounts for the radiative heating properties 
of suspended dust particles and C02 gas. Our goal in developing the model 
is to better understand the processes controlling the exchange of dust and 
water between the surface and atmosphere. 

At last year's DPS meeting, we presented results which demonstrated the 
model's ability to reproduce the basic features of the diurnal variation of 
wind and temperature as measured at the Viking lander 1 site during summer 
(1). We also showed that this diurnal variation was sensitive to dust in 
the atmosphere. Of gr eatest interest in this regard was the diurnal 
variation of the surface stress (Figure 1). It was found that for a given 
wind speed at the top of the boundary layer , the amplitude of the diurnal 
variation in surface stress changed little for modest dust loadings (solar 
optical depths <= 1 ) , but was substantially reduced as the opacity 
approached values typical of mature global dust storms (5 or so). This 
suggests that the change in static stability within the boundary layer due 
to dust heating provides a negative feedback to the dust raising process. 
This, of course, is not a new idea having been suggested earlier by others 
(2,3). However, as far as we know, our calculations represented to the 
first attempt, to quantify the effect. 

A problem with these early calculations is that they were carried out only 
for the Viking 1 lander site during summer, for one assumed value of the 
wind speed at the top the boundary layer (15 m/sec), and for one set of soil 
properties. In addition, the surface heat flux calculation ignored the 
effect of molecular conduction near the surface, an effect that Sutton et 
al . (4) found to be potentially important for Mars. At the Dust III 
workshop, we plan to present results from a variety of simulations that examine 
more fully the behavior of the surface stress and other boundary layer 
parameters, and which include the effects of molecular conduction near the 
surface. The simulations will focus on southern subtropical latitudes 
during the dust storm season. Among the parameters to be varied are dust, 
free atmosphere wind speed, surface roughness, and soil thermal inertia. 
Insight gained from tsese simulations should be useful for interpreting the 
results of the Mars General Circulation Model. 
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FIGURE 1. Model predictions of the diurnal variation of the friction 
velocity at the Viking lander 1 site during early summer for three different 
values of the dust optical depth. The friction velocity is proportional to 
the square root of the surface stress. 
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THE GEOGRAPHICAL VARIABILITY OF MARS SURFACE PHYSICAL 
PROPERTIES. Bruce M. Jakosky, Dept, of Geological Sciences and Laboratory 
for Atmospheric and Space Physics, University of Colorado, and J. Andrew 
Marshall, Colorado Center for Astrodynamics Research, University of 
Colorado. 

Global surface physical properties can be determined by remote- 
sensing techniques such as infrared measurement of thermal inertia and 

radar measurement of fresnel reflectivity. The former is primarily a 

function of the thermal conductivity of the near-surface region, while the 
latter is a function of the dielectric constant and, hence, of the bulk 
density. We have examined the relationship between these two 
parameters using Viking IRTM determinations of thermal inertia and 
Earth-based radar measurements of reflectivity. Data have been binned 
into the same 2° x 2 ° latitude-longitude bins between about -14° and -24° 
latitude, where data are available, excluding data points with significant 

uncertainties. These parameters correlate, although with a correlation 
coefficient of only about 0.5. The general trend of the data suggests that 
conductivity and density are related; several regions are seen to not fall 
along the correlation trend, suggesting other relationships between the 
parameters. 

The general trend can be explained most simply by the spatial 

variation of the degree of formation of a bonded crust in the subsurface. 
An unbonded surface consisting of loose, unconsolidated, fine grains would 
have both a low reflectivity and thermal inertia. Increasing the bonding 

would increase both the conductivity and bulk density of the materials, 
possibly resulting in the observed trend. Alternatively, variations in 
particle size could increase the conductivity while variations in the particle 
shape or degree of packing would produce an increase in density; this 
model is less satisfying in that it requires two properties to vary in 
tandem— size and packing or shape of grains. This latter model cannot be 
excluded, however. 

Several locations do not fall along the general trend; for the most 
part these are geographically well-defined regions which would then have 
distinctive properties. These regions are discussed briefly below. 

Valles Marineris -specifically Capri, Coprates and Eos Chasma, which 
grade into the flood channels to the northeast. These regions have a low 
reflectivity despite having a relatively high thermal inertia. The low 
reflectivity could result either from large wavelength-scale roughness 
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associated with the channel walls and interior or from the presence of 
loose, unconsolidated, but larger-sized particles. The latter interpretation 
is generally consistent with the discussion by Christensen and Kieffer 
(1979) of the relationship between the channel geology and thermal 
properties farther to the north, in the Kasei region. 

Huygens Crater. The walls and interior of the crater have high 
inertia but low radar reflectivity. This is consistent with the idea of 
trapping of moderately coarse material inside the crater, as discussed by 
Christensen (1983). 

Hesparia Planum, to the east of Tyrrhena Patera. This region, 
consisting predominantly, but not exclusively, of ridged plain material, has 
exceptionally high radar reflectivities. Values of 13-18 % suggest bulk 
densities in the range of 1. 6-2.5 g/cm 3 . This value is not typical of ridged 
plain materials globally; values to the west of Tyrrhena Patera have 
reflectivities less than 10 %. 

Cratered terrain. Several locations in the cratered terrain region also 
have reflectivities greater than 13 %, suggesting densities in the range of 
1.75-2.25 g/cm 3 . These points occur as isolated locations and are not 
typical of the cratered terrain units. 

Most regions in the latitude band studied here fit along the general 
trend of thermal inertia with radar reflectivity. Regions having a lower 
reflectivity than expected have been previously identified as unique in 
showing a relationship between the thermal inertia and the underlying 
geology. Presumably, this is related to the deposition and sorting of dust 
and other aeolian debris by the wind, with certain types of surfaces being 
more likely than others to affect the deposition process. Craters and 
channels appear to fit this behavior. Other regions show a higher-than- 

expected reflectivity. These locations are middle-aged geologically rather 
than being younger surfaces; perhaps they have been most-heavily 
weathered during past epochs, with formation of an extensive duricrust, 
and are now being exhumed. 

By examining these regions not fitting the general trend, we appear 
to be identifying the oldest and youngest surfaces at the centimeter scale. 
These surfaces would be interesting ones to sample with a lander or rover 
vehicle. 



POLAR CAP EDGE HEATING IN THE DUSTY NORTHERN WINTER 
ON MARS, Ralph Kahn, Jet Propulsion Laboratory 

There is a qualitative difference between winds 
equatorward of the receding polar cap edge in the southern 
spring, and those in the northern spring. In the south, 
easterly winds begin abruptly before Ls-190°, and continue 
throughout the spring season; in the north, westerly winds 
persist through mid-spring, and easterlies appear only after 
Ls-40° , once the cap edge is at Latitudes poleward of +70° 
(Kahn, 1983). The sense of the near-surface zonal wind at 
these times is set by competition between the meridional 
temperature gradient, which produces westerly acceleration, 
and the CO 2 mass flux from the subliming polar cap, which 
generates easterlies (French and Gierasch, 1979; Haberle et 
al . , 1979). This places an upper bound on the early spring 
CO 2 mass flux from the north cap that is at least a factor of 
50 lower than that produced by the model of James and North 
(1982). 

The symmetric circulation in the northern winter is 
enhanced by the presence of dust, which strengthens this 
circulation (Haberle et al . , 1982) and raises the effective 
atmospheric viscosity (Magalhaes, 1987), moving the high- 
latitude limit of the Hadley cell poleward. Maybe the 
symmetric circulation, and possibly smaller scale waves 
(Barnes and Hollingsworth, 1987) deliver heat that depletes 
the polar cap of CO 2 during the typical northern winter. 
Sudden polar warmings, like the one associated with the 1977b 
global dust storm (Martin and Kieffer, 1979), are unlikely to 
provide the mechanism, since the early spring westerlies are 
found even in years when there is no observed sudden polar 
warming in winter. Spring observations show the cap to be 
patchy equatorward of +70° latitude in several years 
(Christensen and Zurek, 1983; Leovy et al . , 1972 ), and the 
apparent halt in north polar cap regression in early spring 
(James, 1979) could also indicate low frost cover after the 
northern winter. This conjecture suggests a further look at 
high latitude IRTM data for northern winter. 



THE NORTH POLAR SAND SEAS: PRELIMINARY ESTIMATES OF 
SEDIMENT VOLUME; N. Lancaster and R. Greeley, Department of Geology, Arizona State 
University, Tempe, Arizona 85287-1404. 

The North Polar region of Mars contains major accumulations of aeolian sediment in the 
form of extensive sand seas and dunefields (Breed et al., 1979; Tsoar et al., 1979). An estimate 
of the volume of sediment contained in the sand seas is necessary in order to assess their role in 
the martian sedimentary cycle. Spatial variations in sediment volume may provide information on 
the ways in which the sand seas have accumulated, and suggest possible origins and sediment 
sources. In this abstract we provide a preliminary estimate of dune sediment volume based on 
studies of dune morphometry. 

The approach used combines data from terrestrial analogs of martian dunes with detailed 
mapping of dune morphometry (dune spacing, cover) from Viking Orbiter images. Studies of 
dunes in terrestrial sand seas have shown that dune cross-sectional area increases exponentially 
with dune height, whereas dune height increases linearly with dune spacing (Lancaster, 1988). 
Therefore areas of large dunes, although widely spaced, represent a net accumulation of 
sediment. The sediment thickness represented by the dunes is given by their equivalent, or 
spread-out, sediment thickness (EST). The relation between EST and dune spacing for 
terrestrial dunes (Fig.l) provides the means to derive estimates of sediment thickness from a 
parameter (dune spacing) that is relatively easily measured on orbital images. In many areas, 
dunes do not completely cover the surface, so estimates of sediment thickness need to be reduced 
in proportion to the area covered by dunes. These estimates of sediment thickness should be 
considered as minimum values; they are estimates of the sediment contained in the dunes, and do 
not take into account sediment that has accumulated below the dunes. 



Dune spacing (m) 

Fig.l. Relation between EST and dune spacing derived from terrestrial dunes. 

Average dune spacing, and the percentage cover of dunes were measured on Viking 
Orbiter images of the North Polar region of Mars (> 70 0 N.). Calculations of EST were derived 
from these measurements and the relation in Fig.l., and mapped at a scale of 1:5m. (Fig. 3). 
Four major sand seas can be identified (Table 1), similar to those mapped by Tsoar et al. (1979) 
and Dial (1984). Within each sand sea, the area covered by dunes ranges from less than 5% to 
complete cover, with a mean cover of 50.7%. Dune spacing in the sand sea ranges between 150 
and 2955m, with a mean of 512 m (Fig 2). Most North Polar dunes have a crest-to-crest spacing 
of 300-650 m. Very widely spaced dunes (>900 m) are mostly scattered barchans in marginal 
areas of the sand seas. Assuming that martian dunes are a similar shape to terrestrial dunes, most 
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are therefore 12 to 26 m high, based on relations between dune spacing and height for transverse 
dunes in a variety of sand seas (Lancaster, 1988). 





1 


Fig. 2. Distribution of dune spacing (m) in the North Polar Sand Seas 


Equivalent sediment thickness for the sand seas ranges from 0.10-0.50 m to a maximum 
of 5.0-6.0 m, with a mean of 1.81 m (Table 1; Fig. 3). EST is between 2 and 4 m over wide 
areas of the sand sea between 100 and 260° W. Average EST is greatest in this area, and 
significantly lower in the other dune areas. Within each sand sea, EST follows a consistent 
pattern, with the greatest se diment thickness in the center, decreasing toward the margins. 


Table 1. Percentage cover of dunes, mean EST, and sediment volume in major North Polar 
sand seas 


Latitude 

Longitude 

Area 

Dune cover 

Mean EST 

Sediment volume 

(°) 

o 

(km 2 ) 

(%) 

(m) 

(km 3 ) 

74-80 

40-70 

1.06 x 10 5 

37.5 

1.28 

155.6 

76-81 

85-105 

5.00 x 10 4 

42.8 

1.33 

62.3 

76-83 

110-260 

4.70 x 10 5 

57.0 

2.13 

935.7 

76-80 

260-280 

3.25 x 10 4 

33.0 

0.89 

3.3 


Sediment volumes (Table 1) were computed from the isopachs of sediment thickness 
(Fig. 3). The total volume of sediment contained in the dunes of sand seas of the North Polar 
region is 1195 km 3 , of which 78% lies in the large sand sea between 110 and 260° W.. This 
estimate is an order of magnitude less than that of Thomas (1982), but is based upon much more 
realistic views of the size of the majority of martian dunes, and the amount of sand they contain. 

Our new estimate of dune sediment thickness tends to support hypotheses that the dune 
sands were derived from the Polar layered deposits (e.g. Breed et al, 1979; Thomas, 1982). 
There are clear relationships between eroding scarps and areas of dunes. Major areas of dunes, 
such as those between 40 and 70° W, lie downwind of major dissected areas of the perennial 
polar cap. It is perhaps significant that the major area of dunes between 110 and 260° W. lies 
between the present limits and large outliers of the perennial ice. Dissected areas of the polar cap 
and layered deposits extend for some 3.3 x 105 km 2 (ITiomas, 1982). If the layered deposits 
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contained only 1% sand-sized material, some 275 m would have had to have been eroded to 
provide enough sediment to form the dunes. This is consistent with the local relief of 200-800 m 
noted by Blasius and Cutts (1982). 



Fig. 3. Equivalent sediment thickness in the North Polar sand seas of Mars. Isopachs in 
meters. 
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VIKING PHOTOMETRIC AND ALBEDO STUDIES OF 
REGIONAL DUST TRANSPORT ON MARS; S.W. Lee and R.T. Clancy, Laboratory 
for Atmospheric and Space Physics, University of Colorado, Boulder, CO 80309 

There is abundant evidence that, at present, aeolian processes are active over much of the 
surface of Mars (cf. 1; 2). Previous studies have demonstrated that variations in regional albedo 
and wind streak patterns are indicative of sediment transport through a region (3; 4), while thermal 
inertia data [derived from the Viking Infrared Thermal Mapper (IRTM) data set] are indicative of 
the degree of surface mantlir g by dust deposits (5; 6; 7; 8; 9). The visual and thermal data are 
therefore diagnostic of whether net erosion or deposition of dust-storm fallout is taking place 
currently and whether such processes have been active in a region over the long term. 

One of the consequences of variations in regional sedimentation patterns should be relative 
differences in surface texture and macroscopic roughness from region to region; photometric 
analysis of the visual brightness of a surface viewed under a variety of emission and phase angles 
may allow such differences to be quantitatively investigated (10; 11; 12). In addition, 
determination of regional photometric properties will place constraints on the interpretation of the 
related albedo and thermal data sets. 

Disk-resolved (as opposed to disk-integrated) photometric studies of planetary surfaces from 
spacecraft are not common, and photometric studies of Mars using spacecraft data are rare indeed. 
Voyager images of Io have been used to investigate porosity and macroscale roughness of the 
surface (cf. 13; 14); these studies have indicated that it may be more appropriate to study relative 
regional differences in photometric properties, rather than to attempt physically rigorous modelling 
of the results. Images from the early Mariner missions were used to study a few selected regions 
(15); indications were that significant variations in surface properties were detectable. Only the 
first few months of Viking Oi biter IRTM data were systematically analyzed (5); again, significant 
variations in regional properties, as well as the effects of atmospheric aerosols, were detected. 

During the course of the Viking missions, numerous IRTM observing sequences were 
designed with photometric analyses in mind. Several hundred emission-phase-function (EPF) 
sequences, in which the IRTM instrument observed the same area on the surface for an extended 
period as the spacecraft moved overhead, were obtained; of these, about 200 sequences were 
obtained during daytime. Typical sequences consist of several thousand to several tens of 
thousand individual measurements of surface brightness gathered over a wide range of emission 
and phase angles and with spatial resolutions of several tens of kilometers. In addition, about 100 
sequences, scanning the whole disk of the planet from high altitude, were obtained at very low 
phase angles (< 20°). These observations were obtained randomly over much of the planet, but 
most lie within 45° of the equator. The IRTM data set, as a whole, is well calibrated, having been 
corrected for inter-spacecrafi, inter-detector, and temporal calibration variations; a conservative 
estimate for the absolute unc ertainty inherent in this data set is 1-2% (16). Thus, differences in 
regional properties inferred from photometric studies should be attributable to variations in the 
surface and/or atmospheric characteristics and not to calibration errors. To date, however, these 
sequences have not undergone the photometric analysis for which they were intended. 

We are investigating the photometric properties of several regions on Mars through analysis 
of EPF and low-phase-angle IRTM observations. The results of these analyses, in conjunction 
with mapping of regional thennal properties (also from IRTM data) and temporal albedo variations 
(using IRTM and imaging data), are being used to infer details of regional sediment transport and 
consequent properties of the martian surface. The photometric study may also allow quantitative 
investigation of the properties of martian dust and aerosols. In addition to adding to our basic 
knowledge of the surface properties and geologic environment of the martian surface at present, 
these results, particularly those of the photometric analyses, should provide useful input for the 
design of observing sequences by various instruments aboard the Mars Observer spacecraft. 
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CORRELATIONS BETWEEN SURFACE ALBEDO FEATURES AND 
GLOBAL DUSTSTORMS ON MARS: L. J. Martin, Planetary Research Center, Lowell 
Observatory, Flagstaff, A/., and P. B. James, Universiiy of Missouri-St. Louis, St. Louis, 
Mo. 


Seasonal and interarnual variability of albedo features on Mars has been familiar to 
terrestrial astronomers for a long time (Slipher, 1962). The albedo history of the Solis 
Planum area during the Viking Mission has been investigated in detail by S. Lee (1986); 
the period includes the two major 1977 storms observed by the spacecraft. Solis Planum 
was a source for dust during southern spring and summer, when local dust clouds were 
frequently seen in the region, while sedimentation apparently replenished the region 
during the rest of the year. Thus the albedo behavior of dust sensitive regions is prob- 
ably closely correlated with the global dust storm behavior in a particular year. 

The two primary source areas for dust during the growth phases of major dust 
storms on Mars have historically been the Hellas and the Solis Planum regions (Briggs et 
al., 1979). Our 1986 observations showed that both of these areas were unusually dark, 
suggesting a reduced dust inventory in the areas where major storms tend to develop; in 
fact, most of the normally bright features in the southern hemisphere were darker than 
usual in 1986. No major dust storms developed between Lg-160 0 and 1^-330° in 1986, 
a period which includes all of the classic dust storm period on the planet; this supports 
the depletion of dust available in the source regions. During the 1984 opposition, which 
occurred at a somewhat earlier season, Hellas was bright and Solis Lacus was about 
average in size and albedo; however, other albedo features were already darker than 
usual. This suggests that the mechanisms working to deplete these areas, such as the 
local dust storms in the Solis Planum area, were active during 1984 but that there was 
insufficient replenishment from the northern hemisphere to sustain substantial dust activ- 
ity in 1986. 

During most years Hellas can be identified as a large, circular area which is one of 
the brightest features on the planet, comparable to the polar caps and the bright areas, 
such as Arabia, in the northern hemisphere. Hellas was distinctly bright in the dust 
storm years of 1956, 1971, and 1973; but Hellas was no brighter than its surroundings 
on most 1986 images, making its identification difficult. Solis Lacus is normally a dark 
feature surrounded by lighter areas (Thaumasia, Syria, and Sinai); in 1986 Solis Lacus 
was very dark and larger than it had been since 1926. It covered much of the adjacent 
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areas which themselves appeared less bright than in normal years. Most of the southern 
hemisphere, particularly between latitudes -30° and -60°, was darker than usual in 1986; 
the only exceptions were the Ausonia and Eridania regions just to the east of Hellas. 
The behavior of the albedo features in the southern hemisphere during the 1986 appari- 
tion is the exception to the "rule" proposed by Slipher (1962) and others that such areas 
remained light through southern spring and darkened in summer and implies that dust 
inventories were below normal for most of the southern hemisphere in 1986. Thus the 
dust clouds reported in Hellas by some observers (D. Parker, private communication) may 
not have been able to expand into major storms due to a lack of dust. 

If the very dark aspect of the southern hemisphere and the relatively clear 1986 
dust storm season do imply a significant depletion of dust in the source regions, it seems 
likely that mechanisms must exist to return dust from sinks in the northern hemisphere; 
these mechanisms are less well documented than the perihelic duststorms which originate 
in the southern hemisphere. One might speculate that we saw in 1986 an extremum in 
a dust cycle having a period of several Mars years. The 1988 apparition, which spans 
another dust storm season, will give us an opportunity to further investigate the develop- 
ment of these albedo features and their effects on dust storm activity. 
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Lee, S. (1986). LPI Technical Report Number 86-09, 44-45. 

Slipher, E. (1962). Mars : The Photographic Story, Northland Press, Flagstaff, Az. 


45 


EARTHBAS F,D MONITORING OF THE MAFTIAN ATMOSPHERE DUST OPACITY; 
T.Z. Martin, Jet Propulsion Laboratory 

The dust content of the Mars atmosphere is known to exhibit 
considerable interannual variability; specifically, major dust storms 
occur in some years but not others. The effect of dust on geologic 
processes, deposition of polar volatiles, atmosphere dynamics, and remote 
sensing observations makes it important to understand this variability. 
Recent theoretical work on the subject (1) has emphasized the need for 
more data. 

Monitoring of dust behavior on Mars has been sporadic since the failure 
of the Viking 1 lander in 1982. Visual photography once performed 
systematically by the Planetary Patrol has been discontinued. However, 
even when operating they acquired data only when Mars exceeded 13 arc sec 
in diameter - just 20% of the time. Historically, few dust storm seasons 
have been observed from earth or by spacecraft. 

Mars is rarely near opposition during the, prime dust storm season (Ls 
200-300). However, that circumstance obtained in 1986 and 1988. The 
author undertook in 1986 to begin earthbased monitoring of Mars' dust 
opacity in the thermal infrared (2). This work followed naturally the 
analysis performed on the Viking IR Thermal Mapper data to extract 
opacity information. 

The silicate dust band used in the IRTM work (3) falls in the 10 
micrometer window of the earth's atmosphere. Observations in this region 
are unaffected by scattered sunlight, so Mars may be observed in the 
daytime far from opposition. Mars is never less than 4 arc sec in 
diameter, and therefore in principle at least hemispherical differences 
in opacity can be established with mapping performed with a 2 second 
aperture . 

The NASA IR Telescope Facility at Mauna Kea was chosen for its IR 
transparency, excellent tracking and pointing (needed to find stars in 
the daytime), data acquisition software, and especially the facility 
detector system: a SiAs photoconductor coupled with a circular variable 
filter spectrometer covering 6.5-13 micrometers with 1% spectral 
resolution. The silicate band is well resolved without wasting signal, 
and the terrestrial ozone band, also at 9.6 micrometers, can be 
discriminated. A comparison between the appearance of the earthbased 
spectrum and data from the 1969 Mariner 7 IR Spectrometer is shown in 
Fig. 1. 

Derivation of opacity information from the behavior of the silicate band 
rests on the assumption that the dust-bearing atmosphere of Mars during 
the daytime is normally cooler than the surface; the lapse rate is about 
4K/km under clear conditions and about 2K/km when dust content is high. 
In these conditions a silicate absorption appears, whose depth is 
proportional to the dust abundance. While an independent measure of the 
atmospheric temperature would be valuable, such as could be determined 
from the 15 micrometer C02 band, the silicate band behavior alone 
provides a good measure of the changes in opacity. The 15 micrometer band 
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is inaccessible due to the overlapping telluric feature. 

We obtain images of Mars in five wavelengths: 7.8, 8.9, 9.5, 10.4, and 
12.0 micrometers. The first and last of these are "continuum" 
measurements, and the other three fall within the silicate band. 
Conversion to brightness temperature is done using observations of 
standard K type stars lacking silicate shells. From the set of images we 
form a single "band depth" image based on the relative band strength as 
observed in Mariner 9 IRIS spectra. Although absolute calibration of the 
Mars images is a goal, the band depth does not depend strongly on the 
derived absolute temperature. 

Successful measurements have been obtained on ten out of ten observing 
runs completed to date (2-3 half nights each). The seasonal coverage is 
shown in Table 1 . 

At this writing, no NASA funding has been obtained specifically for 
acquisition or analysis of these data; travel support in 1988 was 
provided by The Planetary Society. Reduction of the 750 Mars images 
proceeds slowly. However, detailed analysis of several data sets has been 
done to support proposals. We show in Fig. 2 a single 12.0 micrometer 
image from April 1988, and the corresponding "band depth" image derived 
using data at all 5 wavelengths. The consistency of these derived images, 
which carry model- independent information related to the dust opacity, is 
good evidence that the technique is working. Derivation of dust opacity 
will require assumptions about the emission angle dependence of the 
absorption, the temperature profile, and the surface emissivity. 

The author performed this work as a visiting astronomer at the Infrared 
Telescope Facility, which is operated by the University of Hawaii under 
contract to the National Aeronautics and Space Administration. 
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Table 1 

Seasonal Coverage 


Date 

Ls 

Date 

Ls 

Jul 2-4 

198 

1988 Feb 11-13 

144 

Aug 6-7 

219 

Mar 21-23 

165 

Sep 30-31 

253 

Apr 14-16 

178 

Nov 14-15 

282 

Jun 21-22 

218 



Jul 29-31 

242 

Feb 9-10 

333 

Aug 31-32 

263 

Mar 3-4 

345 

Sep 8-9 

268 



Dec 9-10 

323 


1987 
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Mariner 7 IRS and I RTF Spectra of Mars 
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Fig. 1. Spectra of the silicate band region obtained by the Mariner 7 IR 
spectrometer (top) in 1969 and from the IRTF on June 23, 1988. The IRS 
spectrum has been corrected for instrumental response. The IRTF spectrum 
is ratioed to that of Beta And; a correction for the difference in slope 
between a 300 and 2700K blackbody has also been applied. 
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Fig. 2. Brightness temperature image of Mars, April 17, 1988 1810 UTC, 
and the corresponding silicate band depth image derived from images at 
7.8, 8.9, 9.5, 10.4, and 12.0 micrometers. 



BAND DEPTH, K 







NUMERICAL SIMULATIONS OF GLOBAL DUST STORM DECAY 
J. Murphy, Univ. of Washington, Seattle, WA 98195, O.B. Toon, 

J.B. Pollack, R.M. Haberle, NASA Ames Research Center, Moffett 
Field, CA 94035 

We have been using an aerosol model (1) in one and two dimensions 
to determine as much as possible regarding the physical properties 
of martian dust particles and their suspension lifetimes in the 
martian atmosphere. Our present work involves investigating the 
impact of various physical processes (i.e. sedimentation rate, eddy 
diffusion, coagulation) upon maintaining particles in suspension 
in the martian atmosphere, as well as taking a first look at the 
relationship between atmospheric dynamics and various observations 
of dust in the atmosphere. This work is undertaken in preparation 
for the coupling of the aerosol model and the NASA Ames Mars GCM to 
simulate the atmospheric circulation and dust transport during a 
global dust storm. 

Martian global dust storm decay phases have been observed by 
Mariner 9 (1971) and the Viking landers and orbiters ( 1977a, 1977b) . 
The Mariner 9 IRIS spectra data obtained from the martian southern 
hemisphere subtropics provided information about dust particle 
composition, as well as size and vertical distributions ( 2 ) . These 
data also indicated that the upper and lower portions of the 
atmosphere cleared at the same rate (3), and that the particle size 
distribution did not appreciably change during the observation time 
of (9 100 sols (2). These latter two observations indicate that 
gravitational settling alone was not controlling the rate of dust 
removal from the atmosphere, with vertical eddy diffusion of l.E+07 
cm2/s ( 3 ) and particle coagulation ( 4 ) being suggested as the 
reasons for each of these observations, respectively. 

The visible optical depth values provided by the Viking lander 
imaging experiment (at 23 N latitude), in conjunction with the dust 
particle properties inferred from the Mariner 9 IRIS data, give us 
the opportunity to investigate the importance of the previously 
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mentioned physical processes as they relate to a global dust storm 

decay. With a 1-D version of the aerosol model we have studied 

the impact of these physical processes by reproducing the optical 

depth declines observed at Viking lander 1 during the 1977a & 1977b 

dust storm decay phe ses . The aerosol model contains the previously 

mentioned physically important processes for dust suspension 

maintenance, and by varying their magnitude and/or inclusion in our 

simulations , we arri ved at a combination of these factors which 

produced the best f j t to the data. For initial conditions we chose 

particle size distr j bution # 1 from i 2 ) , constant mass mixing ratio 

with height (3), dust composed of montmorillonite 219b (2), and a 

static, isothermal iitmosphere at 220K. 

The 1977a storm peak optical depth at lander 1 was (§3.2 (5,6), 
decreasing exponentially with a time constant of 75 sols initially, 
and more slowly after 75 sols. In our simulations, spherical dust 
particles fell out much too quickly, even in the presence of strong 
diffusion (l.E+07 cm2/s), to reproduce the observations. Increased 
diffusion and vertical variation of the diffusion had negligible 
additional effect, while the inclusion of Brownian coagulation 
resulted in slightly more rapid optical depth declines. Coagulation 
did little to quench the greater loss of large particles relative to 
the smaller particles. As a means of reducing the Stokes-Cunnungham 
determined particle fall speeds, the particle shape was varied, 
increasing the surface area to volume ratio thereby reducing the 
fall speeds. Disc shaped particles with a thickness to diameter 
ratio of 0.1, in the; presence of vertical diffusion l.E+07 cm2/s, 
resulted in an optical depth decline consistent with the lander 1 
observations of the 1977a storm. Such plate-like particles are not 
inconsistent with those expected for clay particles, and also agree 
with the particle shape resultsfrom lander 1 imaging (5). 

The 1977b storm, with a peak VL1 optical depth of (9 6.0 (6), had 
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an initial exponential decay rate 33% more rapid than the 1977a 

event. This rapid decay rate was best simulated with disc shaped 

particles (thickness/diameter=0. 15) which were a bit "fatter" than 

those of the 1977a storm best fit, and a vertical eddy diffusion 

value of l.E+05 cm2/s. 

This discrepancy between particle shape (and diffusion) for the 
two storms is likely not explained by different particle shapes 
(nor necessarily different diffusion coefficients) for each storm, 
rather we are now investigating a new idea. This idea is that 
dynamics controls the local decay rate of the dust storms by the 
resupply of dust and by vertical velocities that are of the same 
order as the fall velocities. We are in the process of analyzing 
these processes using a 2-D configuration of the aerosol model and 
a 2-D wind field produced by the dynamical model of (7). This 
dynamical model can be run under conditions of varying optical 
depth to produce cross equatorial Hadley type circulations which are 
representative of the atmospheric flow during a global dust storm. 

As well as helping to more clearly define the particle shape 
appropriate for use in 3-D dust storm simulations, we anticipate 
that our 2-D investigations will shed some light upon whether the 
lack of particle size distribution change inferred from Mariner 9 
IRIS data is site specific ( i.e. in the dust source region and/or 
the region of upward vertical motion) or a more global phenomenon. 
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SIMULATIONS OF THE GENERAL CIRCULATION OF THE MARTIAN ATMOSPHERE 
II. DUST STORMS; J.B. Pollack, R.M. Haberle, NASA Ames Research Center, and 
J. Schaeffer, H. Lee, Sterling Software 


We have conducted 
the Martian atmosphere 
equations of meteorolo 
obtain an accurate spe 
dioxide gas at solar a 
scattering calculation 
of the interactions of 
bulk parameterization 
momentum between the s 
Mars consortium inform 
inertia, and albedo of 


numerical simulations of the general circulation of 
with a 3 dimensional model based on the primitive 
gy. Line by line calculations were carried out to 
oif ication of the absorption properties of carbon 
rid thermal wavelengths and single and multiple 
3 were performed to derive an accurate specification 
suspended dust with solar and thermal radiation. A 
scheme was used to evaluate the exchange of heat and 
jrface and the atmosphere. The model incorporated 
ition on the spatially varying topography, thermal 
the surface. 


A large number (16) of numerical experiments were carried out for spatially 
and temporally constant, dust loading (for a given experiment) to determine 
the steady state response of the atmosphere to different choices of dust 
optical depth (0~5) and seasonal date (6 dates spaced about 60° of L^ 
apart). In all cases, a uniform horizontal grid of 7 l 4° of latitude by 9° of 
longitude was used and 13 vertical levels spanning an altitude range of ^7 
km was employed. 


These simulations have a number of implications for dust storm genesis, dust 
transport, and dust removal. Factors that ini luence the magnitude of the 
surface wind stress ard therefore the ability to initiate the saltation of 
sand and the suspension of micron sized dust include seasonal date, the dust 
optical depth, topographical and thermal inertia gradients, and polar cap 
boundaries. Large surface wind stresses occur when Mars is close to 
perihelion, at times cf enhanced dust loading, in association with 
baroclinic wave activity in the winter hemisphere, and in association with 
certain types of topographical slopes at low Latitudes. 

Various components of the circulation respond in different ways to enhanced 
levels of atmospheric dust, thereby influencing both the rate and degree of 
dust spreading from localized sources. For example, both the strength and 
extent of the Hadley circulation are enhanced as the dust loading increases. 


Carbon dioxide condenses in the atmosphere as well as on the surface in the 
winter polar regions. Condensation is likely to take place on suspended 
dust particles and thereby help to remove dust in the winter polar region. 



INFRARED TRANSMISSION MEASUREMENTS OF MARTIAN SOIL 
ANALOGS; Ted L. Roush, NASA Ames Research Center, M/S 245-3, Moffett Field, CA 
94035 

Infrared spectra, in the 2.5 to 25pm wavelength region (4000-400 cm -1 ), of silicate min- 
erals typically exhibit strong absorption features located near 10 - and 20 pm due to the 
stretching and bending modes, respectively, of the SiO 4 tetrahedra ( 1 , 2 ). The complexity 
and wavelength position of these absorptions are used to identify the silicates according to 
structural classes (eg. ortho- versus phyllosilicates) , as well as, specifying its proper position 
within an isomorphous series ( 1 , 2 ). 

Spectroscopic observations of Mars, beyond 5pm which include continuous wavelength cov- 
erage, are limited to data collected by the Mariner 6 , 7, and 9 spacecraft. The first two of 
these missions were fly-bys and covered the ~1. 9-14. 5pm wavelength region (3), while the 
third was an orbiter and covered the 5-50pm wavelength region (4). In each of these data 
sets there is a strong silicate signature near 10pm (5, 6 , 7, 8 ). When Mariner 9 arrived at 
Mars the entire planet was engulfed by a dust storm which slowly cleared with time (7). 
Since the dust particles were suspended, this situation provided an excellent opportunity 
for the comparison of the martian data to transmission measurements of candidate surface 
materials (7). The wavelength position of the 10 pm silicate band has been used to indi- 
cate the presence of material of basic to intermediate composition ( 5 , 6 ), as well as, some 
limited clay compositions (6,7). Figure 1 illustrates brightness temperature spectra from 
the Mariner 9 data and the 10 pm (1000 cm -1 ) band is quite obvious. The Mariner 9 data 
also exhibit an additional strong silicate signature near 20 pm (400-500 cm -1 ) ( 7 ). Previous 
interpretation of the 10 and 20pm bands seen in the Mariner 9 data are most consistent with 
the mineral montmorillonite but some discrepencies remain (7). In particular, the 20pm 
band observed in the martian data exhibit a broad spectral shape which is inconsistent with 
a high degree of silicate polymerization associated with clay minerals ( 7 ), and the relative 
band intensities observed in the martian data ( 10 pm band > 20 pm band) is inconsistent 
with the montmorillonite previously used for comparison purposes ( 9 ). 

Subsequent to the Mariner 9 mission, comparisons of the 0.35-2.5pm earth-based telescopic 
observations of Mars with laboratory reflectance spectra measured for a variety of candi- 
date martian surface materials have resulted in the interpretation that the bright regions 
and soils of Mars are composed chiefly of amorphous (10,11,12) Fe s+ -bearing materials. 

In the visible and near-infrared (0.35-2. 5pm), terrestrial palagonites provide one example 
of spectral analogs for martian surface materials (10,11,12). In the 0.35-2. 0pm wavelength 
region, spectral measurements of iron-substituted smectites provide another possible spec- 
tral analog (13), however beyond 2.0pm their spectra exhibit crystalline clay absorption 
features which are poorly defined in the observed martian spectra ( 12 ). The presence of the 
poorly defined clay bands in the martian spectra has been interpreted as an indication that 
the bright regions and soils on Mars are either slightly more crystalline than the terrestrial 
palagonites, or alternatively, they contain minor amounts of Mg-bearing crystalline clays 
admixed with the palagonites ( 12 ). 

To my knowledge, there have been no previous comparisons of the spectral properties of 
amorphous palagonites at the wavelengths which were covered by the Mariner 9 data. 
Additionally, another useful comparison to the Mariner 9 data would include the spectral 
properties of a variety of smectite minerals. An investigation was undertaken to measure 
the transmission properties of palagonites, as well as, several smectite minerals. While the 
optical constants, necessary for a quantitative comparison to the Mariner 9 data, remain to 
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be derived for these materials, results of the initial transmission measurements are presented 
here and a qualitative comparison to the Mariner 9 data is made. 

Milligram amounts of samples were mixed with KBr powder and then pressed into a coherent 
pellet using approximately 10 tons of pressure. Spectra of the pressed pellets were obtained , 
in the 2.5 to 25/im wavelength region, using a Nicolet 7199 Fourier Transform Spectrometer. 
Each spectrum presented was obtained at a spectral resolution of 4 cm -1 , which remained 
constant at all wavelengths, and consists of 100 individual spectra which were coadded to 
improve signal-to-noise. Spectra presented here were corrected to the same measurement 
conditions but with no sample in the instrument beam. 

Figure 2 shows the spectra of the smectite minerals in the 1800-400 cm -1 region. The 
spectra of these samples all exhibit strong absorptions, near 1000 and 400-500 cm” 1 (10- 
and 20-25 /im), which are due to the silicate fundamentals. A qualtitative comparison to 
Figure 1 illustrates that while in some cases the 10 /im band positions and relative intensities 
of the 10- and 20 ^m spectra of the smectites provide a close comparison to the Mariner 9 
data, the spectra are generally too complex in the 20 /tm region to remain consistent with 
the spacecraft data. 

As seen in Figure 3, the spectra of the two palagonites also exhibit absorptions which 
are centered near 1000 and 400-500 cm” 1 due to the silicate fundamentals. Qualtitative 
comparison of both palagonite spectra to the Mariner 9 data of Figure 1 shows that the 
10 /im band position, the relative intensities of the 10 and 20 /im bands and the relative 
widths of the 20 /tm bands provide a close analog to the Mariner 9 data. However, while 
the spectra of both palagonites provide a better comparison to the Mariner 9 data near 
20 /im than the smectites, they still exhibit more complex structure than that seen in the 
Mariner 9 data. 

In summary, based on silicate absorption band position, relative band intensity, and 20 
/zm band shape and complexity, the 2.5 to 25 /im transmission spectra of two palagonites 
qualitatively provide a better comparison to Mariner 9 data than the spectra of the four 
smectites studied here. However, even the palagonites spectra exhibit a complexity near 20 
/im which is not seen in the Mariner 9 data. 
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Figure 1. Mariner 9 spectra of orbit 56 (top) and 
orbit 8 (bottom), adapted from (7). The gap between 
650 and 800 cm” 1 is due to martian atmospheric CO 2 
absorption. The broad absorptions near 500 and 1000 
cm" 1 correspond to the silicate tetahedral bending 
and stretching fundamentals, respectively. 




Frequency (cm* 1 ) 


Figure 2* Transmission spectra of smectite minerals, 
minimum and maximum values of transmission are in- 
dicated on each curve. Curve A is Ca-Montomorillonite 
(14) (min=57.3, max=84.l). Curve B is nontronite, 
the ferric montmorillonite (15) (mm=40.0, max=87.8). 
Curve C is saponite, the Mg equivalent of montmoril- 
lonite (14) (min=35.0, max=77.7). Curve D is Na- 
montomorillonite (15), (min=16.7, max=81.8). 


Figure 3* Transmission spectra of palagonites, 
weathering products of basaltic glass, VOL15 (15) (top;||||||j! 
min=48.0, max— 84.2) and Pahala Ash (14) (bottom; 
min=48.4, max— 75.5), 
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DUST IN THE POLAR REGIONS OF MARS: IS IT THE SAME AS AT LOW 
LATITUDES? P. C. Thomas, CRSR, Cornell University, Ithaca, NY 

Color and albedo information suggest that some of the northern polar layered deposits are 
similar to the bright areas of Arabia Terra. However, a dark component in the polar deposits 
complicates this interpretation. The mobile component(s) of low latitude variable features is being 
compared to that involved in north and south polar wind streaks. Any significant differences in the 
polar and equatorial dust are important in interpreting the global dust budget. 



DEPOLARIZED RADAR ECHOES AND DUST ON MARS 


T. W. Thompson, Jet Propulsion Laboratory, Pasadena, CA 91109 
H. J. Moore, U.S. Geological Survey, Menlo Park, CA 94025 

We have modeled depolarized echoes of 12.5-cra radio transmissions 
backseat tered from Mars. The model (1) reproduces the variations of the total 
radar cross sections with longitude observed by the Goldstone radar in 1986 
along 7°S (fig. 1); (2) yields larger magnitudes of total radar cross section 
along 22°N than those along 7°S — in accordance with the 1980-1982 Arecibo 
observations (fig. 1); and (3) produces depolarized echo spectra that broadly 
match those observed by Arecibo (fig. 2). In the observations, both polarized 
and depolarized echoes from continuous wave radio transmissions were received 
and sampled as Doppler-spread spectra (1,2, 3, 4, 5). 

In our model, Martian depolarized echoes behave like lunar depolarized 
echoes from the maria, terrae, and young craters. Thus, Martian depolarized 
echo strengths per unit projected area are constant, but echo strengths vary 
from region to region. We divided the surface of Mars into sixty-one regions 
composed of l°x 1° arrays which are called radar map units. These radar map 
units were defined by, first, generalized geologic map units (6,7,8) and, 
second, thermal inertia map units (9,10,11). Depolarized echo strengths were 
then assigned to the radar map units and tested against the observations using 
a "trial and error" method. In the "trial and error" method, the goal was to 
adjust the model echo strengths of the sixty-one radar map units until the 
model (1) matched the magnitudes of the Goldstone radar cross sections at all 
longitudes, (2) produced larger cross sections in the Tharsis region at 22°N 
than at 7°S, and (3) produced depolarized echo spectra that broadly matched 
those observed by Arecibo. 

Model depolarized cross sections of the radar map units vary from 0.007 
to 0.175. Cross sections for cratered uplands and other old terrains are 
0.010, but those for the northern plains range from 0.015 to 0.040. Young 
volcanoes have cross sections that range from 0.070 to 0.175 and those of 
their surrounding lava plains range from 0.050 to 0.160. These cross sections 
are larger than that expected for the average Moon (about 0.005) at the same 
wavelength (1); the cross section of the floor the crater Tycho should be near 
0.07. 


Depolarized echoes are produced by wavelength-size (radii from 0.04 to 
0.4 m) roughness elements at or near the surface and their strengths are 
related to the concentrations of the roughness elements, their depth below the 
surface, and the electrical properties of the matrix containing the buried 
roughness elements (12,13). A semi -quantitative interpretation of the echoes 
(14,5) implies that the concentrations of roughness elements such as slags and 
blocks on lava flows or rocks range from 3% to 76%. 

Large model depolarized echo strengths (0.05 to 0.160) and inferred 
concentrations of roughness elements (22 to 70%) are associated with the lava 
plains of the Tharsis, Amazonis, and Elysium regions where thermal inertias 
are so low that loose dust is required at the surface (9,10,11). Only a few 
centimeters of loose dust are required to produce these low thermal inertias 
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(15) and radar backscatter would be unaffected by such thin layers. On the 
other hand, the low norma L reflectivities of quasi-specular echoes 
(0.041±0.015) from these regions (1,2), experimental data (16,17), and theory 
(18) imply the presence of substantial thicknesses of low density materials 
(about 1,100±200 kg/m 3 ) with relative dielectric constants near 2.310.4. If 
the low density materials, are entirely dusts, local thicknesses of the dusts 
must be at least 3 to 7 m, but only 30 to 78% of the areas need be covered by 
the thick dust deposits. 

Arabia is another region associated with low thermal inertias (9,10,11) 
which again implies that the uppermost few centimeters of materials are loose 
dusts. Here, the model depolarized echo strength (0.008) implies a 
concentration of roughness elements of about 3%; these roughness elements are 
probably rocks because lava flows are rare. In contrast with the Tharsis 
region, the normal reflectivity (0.06110.015) for Arabia (1), experimental 
data (16,17), and theory (18) imply that moderately dense (about 1,400±190 
kg/m 3 ) materials with relative dielectric constants near 2. 7+0. 3 are present 
within the uppermost 2 to 3 m. Because the radar waves would be unaffected by 
a few centimeters of loose dust, the surface materials of Arabia may be 
chiefly (97%) moderately dense materials overlain by a few centimeters of dust 
and rosks with radii from 0.04 to 0.4 m (3%). The moderately dense materials 
may be cemented or compacted dusts. Thus, the rate that processes cement and 
compact dust in Arabia may, in fact, be comparable to the rate that loose dust 

accumulates (19). 

Analyses of both the polarized and depolarized echoes from continuous 
wave radio transmissions at 12.5-cm wavelength are in progress and additional 
observations will be made this year. These observations and analyses will 
provide important information on the distributions, thicknesses, and nature of 

dusts on Mars. 
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DEPOLARIZED RADAR ECHOES 
Thompson, T. W. and Moore, H # J 



O ARECIBO 1900-1982 OBSERVATIONS □ MODEL: ARECIBO 1980-1982 
Him GOLDSTONE 1986 OBSERVATIONS »-»MODEL: GOLDSTONE 1986 


Figure 1. Model and observed total radar cross sections for depolarized 
echoes as a function of longitude. Goldstone cross sections are averages in 
10° longitude bins between 3° and 14° S. Arecibo cross sections are single 
observations between 22° and 25° N. Model cross sections are at 7°S and the 
same latitudes as the observations in the north. 
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Figure 2. Model and observed depolarized echo spectra for two Arecibo 
observations at the indicated latitudes and longitudes. Note that the model 
and observed spectra are similar. The double hatched area shows the model 
spectrum for a featureless, uniformly bright planet with a total cross section 
of 0.010. K and L represent spectra dominated by lava plains and volcanoes. 

0 and P represent spectra dominated by cratered uplands and Arabia. 





MARS GLOBAL ATMOSPHERIC OSCILLATIONS: ANNUALLY SYNCHRONIZED, 
TRANSIENT NORMAL MODE OSCILLATIONS AND THE TRIGGERING OF GLOBAL DUST 
STORMS; J. E. Tillman, Department of Atmospheric Sciences, AK-40, University of Washington, 
Seattle, Wa., 98195 

Observation and analyses of the atmospheric pressure at the surface on Mars for several 
Martian years, has revealed transient events in the daily pressure variations which may be Kelvin 
normal mode oscillations and these or similar oscillations may be a necessary ingredient in the 
initiation of Martian global dust storms. The transient events last only a few sols, ( Mars days ), 
appear to repeat on an annual basis to within a few sols, appear to cover a large portion of the 
hemisphere on the same sol or at least within a few sols, occur in pairs with separations of 20 sols 
( in some instances ), and are coincident with the annual pressure minimum, ( i.e., the maximum 
concentration of carbon dioxide in the southern hemisphere ). Prior analyses of these transient 
pressure oscillations suggested that they differed from exact harmonics of the diurnal cycle but 
did not have adequate spectral resolution to make any definitive conclusion. A reanalysis of the 
data with higher resolution shows that they consist of spectral components generally near to and 
sometimes identical in frequency with the diurnal and semi-diurnal harmonics. It is suggested that 
these events are transient, Kelvin, normal mode, global oscillations due to their close agreement in 
frequency with theoretical calculations. Simulation of composite diurnal and almost diurnal har- 
monics and analyses with the same techniques and parameters demonstrates that the two nearby, 
but distinct, spectral lines apparently are not artifacts of the analytical techniques. Comparison 
of the frequencies with theoretical calculations of Hamilton and Garcia [1], 1984 and Zurek, (pri- 
vate communication 1988), suggests that these are Kelvin, wavenumber 1 and 2 modes which are 
symmetric with respect to the equator. 

Examination of the initiation of the 1977 A, 1977 B and 1982 A global dust storms indicate that 
pseudo-diurnal harmonic oscillations might be involved in triggering the global dust storms, possibly 
in conjunction with winds due to baroclinic and Hadley circulations. Prior to the 1977 B global dust 
storm, the diurnal harmonic disappears and is replaced by an oscillation at 1.1 cycles/sol. If the 
lack of a classical diurnal tide indicates the suppression of convection, this might partially explain 
the slow growth of the 1977 B storm. The daily average pressure during the transient season, which 
preceeds the global dust storm season by almost 1/4 of a year, exhibits a different behavior during 
the global dust storm years than during non global dust storm years. This implies that the global 
circulation must differ significantly in these cases and may be related to the presence or absence of 
global dust storms later in a given year. 

Wind data will be discusse d in the context of changes during the transient events to assess their 
support for the Kelvin mode hypothesis. In general, these analyses suggest that normal mode, high 
frequency global oscillations, different from the classical solar driven tides, may be quite common 
on Mars and may be necessary for the initiation of the global dust storms. 
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MARTIAN DUST THRESHOLD MEASUREMENTS - SIMULATIONS UNDER 
HEATED SURFACE CONDITIONS 

White, Bruce R., Dept, of Mechanical Engineering, University of California, Davis, CA, 
95616, and Ronald Greeley, Dept, of Geology, Arizona State University, Tempe, AZ 85287. 


On Mars, diurnal changes in solar radiation set up a cycle of cooling and heating of the 
planetary boundary layer which strongly influences the wind field. Early in the morning, before 
sunrise, the stratification of the air layer is stable because the ground is cooler than the air above it. 
When the sun rises, solar radiation heats the ground faster than the air above it. When the ground 
is warmer than the air, heat is transferred to the air above the surface and parcels of the heated air 
move upward past the cooler air. On clear days the air receives almost all its heat from the ground. 
The heated parcels of air near the surface will, in effect, increase the near surface wind speed or 
increase the aeolian surface stress the wind has upon the surface, when compared to an unheated or 
cooled surface. This means that for the same wind speed at a fixed height above the surface, 
ground-level shear stress will be greater for the heated surface (unstable case) than an unheated 
surface. Under these conditions, it is possible to obtain saltation threshold at lower mean wind 
speeds than when the surface is not heated. 

The boundary layer is in neutral stability when the wind profile is exactly logarithmic. The 
simplest interpretation of the eddy structure would be a set of circular eddies with diameters given 
by the mixing length 1 = kz, which rotate at a uniform tangential speed equal to the friction velocity 
u*; i.e., w' = u' = u* = 1 9u/9z, where w' and u' represent the respective vertical and horizontal 
turbulent velocity fluctuations in the flow. In unstable conditions, vertical motions are enhanced 
by buoyancy due to an extent inversely dependent on the existing wind shear 9u/9z. In terms of 
velocity fluctuations, w’ exceeds u', i.e., 

w' = u' + [a contribution generated by buoyancy] 
where w’ is given by 1 9u/9z, but 1 is greater than kz. The converse occurs in stable (inversion) 
conditions. The Richardson number is proportional to the ratio of flow energy derived from 
buoyant air forces to flow energy obtained from the shear of the wind over the surface. The 
Richardson number provides quantitative data on the destabilizing effect of buoyancy and the 
stabilizing effect of large-scale wind. The Richardson number is 



The Richardson Number may be approximated by the bulk Richardson Number which is given by 

_ , g A0Z 

Ri = -* *- 


where A0 is the potential temperature difference between the ground and height Z, T is the mean 
temperature between the ground and height Z, and U is the mean wind speed at height Z. 

In order to assess the effect of surface heating on Mars and its influence on particle 
entrainment, threshold measurements have been made in MARSWIT. This is an open circuit wind 
tunnel that is operated within a chamber that allows a range in operating pressure from a few 
millibars (appropriate for Mars) to one bar of pressure (Greeley et al., 1977). Experiments were 
carried out with both heated and unheated surfaces. The heated surface represents diurnal surface 
heating by radiation from the sun. The unheated surface represents a neutrally stable condition. 

The heating of the floor will model solar radiation that affects the vertical turbulent structure of the 
boundary layer flow. Although exact heating is unknown for Mars, it is expected to produce a 
maximum surface - freestream temperature difference of 25 K (Hess et al., 1977; Ryan and Henry, 
1979). Dust (~11 pm) was placed on the tunnel floor by aerodynamic settling (to model the natural 
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emplacement of dust from the atmosphere). The test bed was approximately 40 cm by 100 cm. 
Threshold was determined by observation via closed circuit television and by an electrometer 
particle detector. Vertical temperature profiles were measured under several heating conditions to 
determine friction speed, u* as a function of ffeestream speed at a specified heating condition. 
Additionally, the surface dust temperature was measured with a thermocouple from which a value 

of bulk Richardson number was determined. . 

Although there are no direct data available as to the mean wind speed and surface temperature 
when dust movement occurred on Mars, Arvidson et al. (1983) estimate that winds of 25 to 30 m/s 
would be needed to initiated particle motion in loose soil. Presumably this would occur at Mars 
noon when the temperature difference between the surface and Lander was a maximum, and 
corresponds to a bulk Richardson number of about -0.002 as a threshold value. Wind tunnel 
results show particle movement occurs in the bulk Richardson number range of -0.0012 to 
-0.0023, when scaled to the height of the Viking Lander (1.61 m), well within the value estimated 
to occur on Mars during the initiation of dust movement (personnel communication, R. Haberle). 
These wind speed values reduce the mean wind speed needed to initiate dust motion by 25 to 40%, 
depending upon the exact value of surface heating or value of bulk Richardson number, as 
compared to an unheated surface. 
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MARTIAN GREAT DUST STORMS: APERIODIC PHENOMENA? R. W. Zurek, 

Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 
91109, and R. M. Haberle, NASA Ames Research Center, Moffett Field, CA 94035 

The frequency of occurrence of martian great dust storms accounts for 
much of the interannual variability of the atmospheric thermal structure 
and circulation on Mars. These storms profoundly affect the seasonal 
climatic cycles in any year in which they occur. On very long (i.e., 
astronomical) time scales, the great dust storms may provide the dust which 
appears in the polar laminated terrain and in long-term dust deposits 
elsewhere on the planet. 

On all these time scales, from seasonal to interannual to very-long- 
term, the frequency with which dust storms occur is critical to quantifying 
their climatic effects. There are two limiting scenarios for this frequency 
of great dust storm occurrence: (1) the great dust storms occur nearly 

every martian year, in response to a seasonally occurring "trigger," or (2) 
the great dust storms are not periodic, or even quasi-periodic, at all, but 
rather occur as part of a highly nonlinear, stochastically forced system. 

Zurek (1982) argued that, based on the Viking Lander surface pressure 
data and on both Earth-based and spacecraft data since 1956, a great dust 
storm (defined as one which is of planetary scale, though it certainly need 
not be planet-wide or even planet-encircling) could have occurred nearly 
every Mars year. Furthermore, it was argued that this scenario was 
consistent with Earth-based observations earlier in this century, given the 
natural limitations of viewing Mars from Earth. 

Martin (1984) countered by pointing out that relatively few planet- 
encircling dust storms had ever been viewed on Mars and that these had all 
occurred since 1955. In particular, he argued that the absence of great 
dust storms during earlier periods was due to a real lack of storms and not 
to viewing limitations. 

Leovy et al. (1986) noted the two very different northern winter 
weather regimes that occurred even during the Viking years, depending on 
whether or not a planet-encircling dust storm occurred that year. They 
speculated that the occurrence of the storms themselves might be an example 
of a bifurcation phenomenon whose alternate paths are determined by randomly 
occurring natural variations in the atmosphere-dust-polar cap system during 
earlier seasons. On the other hand, Haberle (1986) suggested that the 
occurrence of the great dust storms was somewhat deterministic, in that the 
occurrence of a great dust storm one year might suppress the generation of 
such a storm the following year. 

Dust accumulation in surface sinks over time would be very different 
under these different scenarios for the frequency of the great storms. The 
present "dusty period" of frequent great dust storms may be quite anomalous. 
If so, the amount of dust carried into the polar regions during long-term 
cycles could be much less than that estimated previously from the annual 
dust cycles observed in the last few decades (e.g., Pollack et al., 1979). 
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After a brief review of current hypotheses for dust-storm generation 
and the frequency of that generation, we investigate a model for dust-storm 
occurrence that would yield a highly aperiodic time series. Following 
Leovy et al. (1986), we envision that there are two quasi-stable climatic 
states for Mars, corresponding to years with and without great dust storms, 
and that the Mars climate system jumps between these two states due to 
randomly occurring natural perturbations. However, we hypothesize that 
this transition is asymmetric — that is, Mars may spend more time in the 
non-dust-storm regime, than in the dusty regime — because of nonlinearities 
in the physical system. For Mars, the nonlinear feedback may be provided 
by the radiative-dynamic effects of airborne dust hazes. 

The specific model is motivated by a terrestrial model (Vallis, 

1986) in which El Nirfo-Southern Oscillation climatic fluctuations are 
simulated as a chaotic dynamical system. These terrestrial phenomena 
resemble the great dust storms on Mars to the degree that both may be 
aperiodic episodes which have a preferred season for onset and that the 
strength of the response (as measured by sea-surface temperatures on Earth 
or visible opacity on Mars) may be quite varied, as will be the global 
climatic effects. 
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MARS OBSERVER PMIRR: DISTRIBUTION AND 

TRANSPORT OF ATMOSPHERIC DUST; R. W. Zurek, R. D. Haskins, J. 
T. Schofield and D. J. McCleese, Jet Propulsion Laboratory, 
California Institute of Technology, Pasadena, CA 91109. 

The Pressure Modulator Infrared Radiometer (PMIRR) is a 
nine-channel visible and thermal infrared radiometer designed 
to make the measurements needed to retrieve vertical profiles 
of atmospheric temperature, airborne dust extinction, and 
atmospheric water vapor, and to distinquish between CO 2 and 
H 2 O condensate hazes and clouds, all as a function of 
atmospheric pressure. Data from PMIRR's single broad-band 
(0.3 - 3.0 pm) solar channel and eight thermal infrared 

channels will also be used for radiation budget studies. 
Geostrophic wind shears are constructed from the observed 
temperature fields, while zonal-mean winds are estimated from 
the zonally averaged thermodynamic energy equation, using 
observed temperatures and the diabatic forcing calculation 
from airborne dust observations. 

With regard to dust, the PMIRR measurements will provide 
a three-dimensional global picture of the t ime - e vo Iving 
atmospheric dust loading and a zero-order estimate of the 
role of atmospheric transport in that evolution. Pre-launch 
activities to be carried out by the PMIRR Science Team and 
collaborators include the development and testing of fully 
integrated atmospheric retrieval algorithms and the 
assessment of different methods for using the PMIRR 
measurements to elucidate the atmospheric transport of dust. 
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